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was  developed,  new  display  methods  using  color  were  Implemented  to  present 
the  data,  F-reglon  analysis  procedures  were  extended,  and  binning  procedures 
were  developed  to  facilitate  statistical  analysis  of  the  data.  In  addition, 
workshops  were  held  so  that  the  observations  could  be  discussed  and  inter¬ 
preted  by  the  MITHKAS  participants. R  As  the  studies  developed,  the  results 
have  been  presented  at  meetings  and  k  large  number  of  publications  have 
resulted.  Because  it  was  such  a  valuable  data  set  and  as  a  result  of  the 
above  activities,  the  number  of  MITHRAS  participants  grew  considerably  from 
those  actually  supported  by  AFOSR.  Although  many  studies  have  been  com¬ 
pleted,  others  begun  later  or  based  on  earlier  results  are  under  way. 

MITHRAS  has  become  the  model  for  other  global  studies  involving  the  full 
network,  of  six  radars.  \ 

The  project  has  led  to  the  following  scientific  accomplishments: 

The  best  view,  to  date,  has  been  obtained  of  the  instantaneous  shape  of  the 
convection  pattern  in  the  auroral  region  and  polar  cap  by  combining 
Millstone  Hill,  Chatanlka,  and  DE-2  electric-field  measurements.  A  series 
of  case  studies  established  that  the  convection  intensifies  simultaneously 
in  the  morning  and  afternoon  convection  cells.  During  a  substorm,  another 
case  study  showed  that  the  large-scale  convection  pattern  could  remain 
stable  for  a  prolonged  period,  even  though  the  period  was  one  of  sustained 
magnetic  activity.  Average  patterns  of  plasma  convection  were  calculated 
for  two  different  orientations  of  the  IMF  By  component.  The  patterns  ob¬ 
tained  showed  that  the  effect  of  the  By  component  is  stronger  and  different 
than  previously  thought.  At  or  near  the  convection  reversal  in  the  after¬ 
noon  cell,  l.e.,  at  the  polar  cap  boundary,  there  can  be  large  energy  inputs 
to  the  electron  gas  by  soft  particles  and  a  large  downward  heat  flux.  The 
winter  nighttime  F-region  density  strongly  depends  on  longitude;  this  de¬ 
pendency  was  attributed  to  cross-polar  cap  transport  of  ionization  and  to 
the  fact  that,  because  of  the  offset  of  the  geographic  and  geomagnetic 
poles,  the  solar  illumination  of  the  large-scale  convection  pattern  varies 
with  UT .  Using  more  inputs  and  being  subject  to  more  constraints  than  ever 
before,  a  theoretical  model  of  the  blgh-latltude  ionosphere  was  able  to 
successfully  reproduce  auroral-zone  observations  made  from  two  widely  spaced 
longitudes.  Theoretical  calculations  of  the  6300-A  intensity  from  atomic 
oxygen,  showed  that  the  hot  electrons  found  near  the  afternoon  convection 
reversal  are  responsible  for  a  much  larger  proportion  of  the  emission  than 
previously  considered.  Although  convection  has  the  dominant  role  in  deter¬ 
mining  the  large-scale  nighttime  meridional  neutral  wind  at  F-reglon 
altitudes,  local  effects  such  as  particle  and  Joule  heating  can  create 
significant  perturbations.  At  midnight  in  the  polar  cap  a  new  phenomenon 
was  observed,  an  abatement  or  sometimes  even  a  reversal  of  the  southward 
meridional  wind.  The  importance  of  ion  drag  was  further  demonstrated  in 
statistical  analyses  that  showed  the  meridional  component  of  the  neutral 
wind  to  be  significantly  correlated  with  the  By  component  of  the  IMF.  These 
analyses  also  showed  that  the  neutral  wind  varies  significantly  with  season. 
Several  phenomena  observed  in  the  early  morning  at  Chatanlka — such  as 
elevated  ion  and  neutral  temperatures,  significantly  perturbed  meridional 
neutral  winds,  and  reduced  plasma  densities — could  be  the  result  of  very 
large  Joule  heating  Just  poleward  of  the  radar  or  in  the  polar  cap. 

Finally,  new  methods  were  devised  to  determine  the  neutral  exospheric  tem¬ 
perature  from  high- latitude  radar  data;  the  daily  and  solar-cycle  variations 
were  studied. 
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I  INTRODUCTION 


A.  The  MagnetospherlcrIanosBlierlc-Thermg.s.Pherlg_SY.stgm 

Since  the  realization  in  1958  that  the  solar  wind  exists,  a  model 
has  developed  of  the  outer  limits  of  our  environment  and  the  interact¬ 
ions  among  its  components.  The  solar  wind,  dragging  the  solar  and 
interplanetary  magnetic  field  (IMF)  with  it,  interacts  with  the  earth's 
magnetic  field  to  form  the  comet-shaped  magnetosphere,  Figure  1-1.  It 
also  sets  up  a  dawn-to-dusk  electric  field  across  the  tail  of  the  magne¬ 
tosphere,  and  it  helps  to  populate  the  magnetosphere  with  charged  part¬ 
icles.  The  magnetosphere  is  closely  associated  with  the  ionosphere  and 
neutral  atmosphere  at  high  latitudes  because  of  the  topology  of  the 
magnetic  field.  The  polar-cap  field  lines  extend  to  the  magnetopause 
(the  outer  boundary  layer)  and  out  into  the  tail  lobes.  The  auroral- 
region  field  lines  extend  to  the  boundary  plasma  sheet  and  central 
plasma  sheet  on  the  nightside  and  to  the  cleft  on  the  dayside.  The 
field  lines  Just  equatorward  of  the  auroral  region  connect  to  the  ring 
current.  These  associations  lead  to  strong  interactions  among  the 
magnetosphere,  ionosphere,  and  neutral  atmosphere.  Most  of  which  are 
concentrated  in  the  D,  E,  and  F  regions  of  the  ionosphere,  and  the 
mesospheric  and  thermospheric  portions  of  the  neutral  atmosphere. 

Because,  to  a  first  approximation,  the  magnetic  field  lines  are 
equipotentials,  the  electric  field  in  the  magnetosphere  maps  along  the 
field  lines  into  the  ionosphere,  causing  the  ambient  electrons  and  ions 
in  the  northern  and  the  southern  high-latitude  regions  to  drift  such 
that  two  convection  cells  are  formed.  Figure  1-2,  that  cover  the  polar 
cap  and  auroral  region.  The  ions,  thus  set  into  motion,  collide  with 
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FIGURE  1-1  SCHEMATIC  REPRESENTATION  OF  THE  MAGNETOSPHERE 


(Source:  Heppner,  19771 


FIGURE  I-: 


IF 

EMPIRICAL  MODEL  OF  THE  IONOSPHERIC  CONVECTION  CELLS  IN  THE 
NORTHERN  HEMISPHERE.  They  are  shown  in  magnetic  coordinates  for  an 
observer  on  the  earth.  The  contours  are  for  electrostatic  potential  with  4-keV 
steps.  The  dashed  line  is  the  Harang  discontinuity  . 


neutrals,  heating  both  the  neutrals  and  Ions  (Joule  heating)  and 
transferring  momentum  to  the  neutrals  (ion  drag).  The  neutrals  are  set 
into  motion  as  a  result  of  both  the  Joule  heating  and  the  ion  drag.  In 
the  E  region,  the  ion  motion  is  impeded  by  these  ion-neutral  collisions, 
while  the  electrons  suffer  very  few  collisions.  The  resultant  dif¬ 
ference  between  ion  and  electron  velocities  is  the  origin  of  the  auroral 
electro jet  currents,  to  which  magnetometers  respond.  The  same  velocity 
difference  further  leads  to  the  development  of  instabilities  in  the 
plasma,  which  give  rise  to  the  density  irregularities  used  for  the  STARE 
and  SABRE  radars.  In  addition,  the  instabilities  give  rise  to  large 
electron  heating.  Energetic  particles  (mostly  electrons)  following  the 
field  lines  from  the  cleft,  central  plasma  sheet,  and  boundary  plasma 
sheet  precipitate  into  the  atmosphere,  forming  the  auroral  oval, 

Figure  1-3,  in  each  high-latitude  region.  These  precipitating  particles 
deposit  energy  in  the  atmosphere  at  D-,  E-,  and  F-region  altitudes.  In 
addition  to  exciting  neutrals  and  ions,  producing  the  auroras,  the 
precipitating  particles  ionize  the  neutral  gas,  giving  rise  to  electron 
density  enhancements,  as  well  as  heat  the  ambient  electrons  and  neu¬ 
trals.  The  latter,  like  Joule  heating,  affects  the  neutral  wind.  In 
addition  to  the  energetic  particles,  thermal  particles  flow  in  both 
directions  along  the  field  lines.  These  energetic  and  thermal  particles 
carry  field-aligned  (Blrkeland)  currents.  Figure  1-4,  that  further 
couple  the  ionosphere  and  magnetosphere.  These  currents  close  in  the 
ionosphere  and  in  the  magnetosphere. 

Many  of  these  interrelations  are  presented  schematically  in  the 
auroral  diagram.  Figure  1-5.  The  geophysical  parameters  are  shown  in 
circles  and  rectangles  with  rounded  corners.  The  physical  processes 
that  affect  them  are  shown  in  rectangles  (with  square  corners). 

Transfer  of  energy  and  momentum  from  the  magnetosphere  occurs 
continuously  and  periods  exist,  called  substorms,  during  which  the 
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[Source:  I ijima  and  Potemra,  1978) 


FIGURE  1-4  SCHEMATIC  REPRESENTATION  OF  THE  LARGE-SCALE  BIRKELAND 
CURRENTS 


transfer  increases  impulsively  and  remains  high  for  several  hours.  That 
the  energy  transfer  alone  can  have  major  effects  on  the  ionosphere  and 
neutral  atmosphere  is  apparent  when  we  realize  that  the  energy  deposited 
above  90  km  in  a  localized  area  can  exceed  that  which  would  be  deposited 
by  the  overhead  sun  by  a  factor  of  10  to  100.  Indeed,  large  effects  do 
occur  in  every  geophysical  parameter:  electron  densities,  electron  and 
ion  temperatures,  neutral  winds,  neutral  temperatures  and  densities. 

In  addition  to  these  auroral  phenomena  in  the  high-latitude  upper 
atmosphere,  this  part  of  the  atmosphere  already  behaves  very  differently 
from  midlatitudes  because  much  of  it  is  above  the  arctic  circle  where 
it  can  have  as  much  as  24  hours  of  darkness  or  sunlight,  depending  on 
the  season.  This  variable  insolation,  for  instance,  leads  to  major 
seasonal  and  hemispheric  differences  in  the  F-region  densities  and  the 
E-region  conductivities.  Thus  one  of  the  challenges  in  this  region  is 
to  distinguish  between  effects  that  are  driven  by  or  modulated  by  the 
solar  wind  and  those  that  are  driven  or  modulated  by  the  solar  EUV  and 
UV  radiation. 

The  interactions  among  the  magnetosphere,  ionosphere,  and  thermo¬ 
sphere  lead  to  a  vast  array  of  phenomena,  which  are  of  great  scientific 
interest.  In  addition,  the  changes  that  result  in  the  ionosphere  and 
thermosphere  can  have  major  effects  on  a  variety  of  systems  that  involve 
the  polar  regions. 

During  periods  of  increased  geomagnetic  activity  the  best  known 
phenomenon  is  the  aurora.  The  emissions  occur  not  only  in  the  visible 
but  also  in  the  ultraviolet  and  infrared  parts  of  the  spectrum.  They 
can  present  a  confusing  background  that  hampers  locating  other  emission 
sources.  Almost  as  well  known  are  many  radio  or  propagation  phenomena. 
These  include  the  disruption  of  HF  radio  communications,  which  also 
implies  difficulty  for  direction  finding  and  over-the-horizon  radar.  At 


higher  frequencies,  these  phenomena  include  the  fading  and  scintillation 
of  satellite  signals.  They  also  include  variations  in  total  electron 
content,  on  various  time  scales,  that  affect  satellite  ranging.  The 
neutral  density  at  high  altitudes  is  increased,  giving  rise  to  increased 
drag  on  satellites.  In  the  mesosphere,  ozone  is  reduced.  Speculation 
exists  that  tropospheric  weather  can  be  affected.  Currents  are  induced 
in  the  lower  ionosphere  and  at  ground  level.  The  former  gives  rise  to 
magnetic-field  fluctuations,  the  latter  to  potentially  destructive 
currents  on  communication  lines  and  power  grids. 

B.  MITHRAS 

The  purpose  of  Project  MITHRAS  (an  acronym  for  Magnetosphere  Iono¬ 
sphere  THermosphere  RAdar  Studies)  is  to  increase  our  understanding  of 
these  Interactions  and  the  resulting  effects.  The  observations  were 
made  during  a  brief  period  when  three  radars  could  observe  much  of  the 
auroral  region,  i.e.,  much  of  the  auroral  oval  and  equatorward  portion 
of  the  convection  pattern. 

The  individual  radars  are  each  very  powerful  research  tools.  An 
extensive  set  of  ionospheric  and  neutral  atmospheric  parameters  can  be 
determined  from  the  radar  measurements.  Many  of  these  parameters  are 
mentioned  in  Section  III,  in  which  we  discuss  the  research  results. 
Numerous  altitudes  can  be  probed,  yielding  profiles  of  these  parameters. 
The  fact  that  altitude  profiles  are  obtained  is  important  because  what 
is  observed  at  one  altitude  is  often  related  by  transport,  diffusion,  or 
conduction  to  what  happens  at  other  altitudes.  Moreover,  the  altitude 
range  observed  at  most  radars  can  include  the  region  between  80  and 
150  fcm  that  is  inaccessible  to  satellites  and  only  accessible  to  rockets 
for  a  few  minutes.  The  high-latitude  radars  can  probe  a  wide  spatial 
area  within  seconds  or  minutes.  Thus,  for  instance,  they  can  establish 
the  location  of  a  measurement  with  respect  to  a  convection  or  precipita- 


tion  boundary.  Finally,  they  can  determine  the  time  evolution  of  the 
parameters.  Individual  profiles  require  typically  5  s  to  1  min.  How¬ 
ever,  an  experiment  usually  consists  of  many  profiles  in  different 
directions,  leading  to  an  experiment  time  resolution  of  5  min  to  40  min. 
A  radar  can  probe  the  same  region  for  periods  that  typically  range  from 
8  hrs  to  72  hrs.  In  contrast,  a  satellite  will  return  to  the  same 
latitude  about  every  90  min.  Depending  on  the  orbit,  a  satellite  may 
take  weeks  to  sample  24  hrs  of  local  time  (thereby  mixing  together  many 
effects  such  as  season  and  geomagnetic  activity),  or  it  may  stay  fixed 
in  local  time  (LT). 

Despite  the  great  capability  of  an  individual  radar,  a  number  of 
problems  exist  that  can  best  be  addressed  simultaneous  observations  with 
by  as  many  radars  as  possible.  When  one  radar  detects  a  change  in  a 
parameter,  such  as  convection  velocity,  as  a  function  of  time,  it  is  not 
clear  whether  the  change  arose  because  the  phenomenon  evolved  [a 
variation  in  universal  time  (UT)]  or  because  the  radar  moved  with 
respect  to  a  steady-state  pattern  with  spatial  gradients  [a  variation  in 
local  time  (LT)].  Because  LT  is  equivalent  to  longitude,  this  quandry 
is  also  referred  to  as  a  space-time  ambiguity. 

Another  situation  that  introduces  LT-UT  differences  is  that  the 
high-latitude  phenomena  are  affected  by  two  sources  of  energy  input: 
the  solar  UV  and  EUV  energy  input  depend  on  local  time  and  geographic 
latitude,  or  equivalently,  longitude  and  latitude;  the  magnetospheric 
energy  and  momentum  transfer  depend  on  local  magnetic  time  and  geomag¬ 
netic  (invariant)  latitude.  Because  the  geographic  and  geomagnetic 
poles  do  not  coincide,  these  coordinate  systems  do  not  coincide.  The 
result  is  that  at  each  location  on  the  earth,  there  is  a  different 
mixture  of  inputs  related  to  these  two  energy  sources.  Observationally, 
this  means  that  sites  at  the  same  latitude  (geographic  or  geomagnetic) 
will  measure  different  values  corresponding  to  different  mixtures  of 


inputs.  Accordingly,  it  also  means  that  a  UT  variation  will  be  imposed 
on  observations  made  at  a  given  site  as  the  geomagnetic  pole,  hence 
auroral  oval,  moves  about  the  geographic  pole. 

The  first  phase  of  MITHRAS  Cde  la  Beaulardlere  et  al..  1982,  1984b] 
was  an  observation  program,  between  May  1981  and  June  1982,  based  on 
three  incoherent-scatter  radars  that  could  observe  the  northern  high- 
latitude  region  between  about  57°  and  73°  invariant  latitude  (A): 
Chatanika,  Millstone  Hill,  and  EISCAT  (the  European  Incoherent-SCATter 
radar).  These  latitudes  typically  embrace  the  auroral  oval  in  the 
midnight  sector  and  are  equator ward  of  it  in  the  noon  sector.  This 
period  was  unique  because  three  such  radars  were  at  almost  equally 
separated  longitudes.  Figure  1-6.  It  began  Just  before  the  beginning  of 
observations  at  EISCAT  and  ended  shortly  after  the  dismantling  of  the 
Chatanika  radar  for  its  move  to  Sondre  Stromfjord,  Greenland.  During 
this  limited  period,  the  AFOSR  support  enabled  roughly  three  simul¬ 
taneous  24-hr  observations  per  month  for  10  months.  (The  catalog  of 
observations  is  given  in  Appendix  A.)  As  seen  in  Figure  1-6,  the  new 
configuration  of  Millstone  Hill,  Sondrestrom,  and  EISCAT,  starting  in 
1983,  is  very  different.  With  the  radars  closer  together,  it  is  more 
suited  for  detailed  studies.  It  also  differs  in  that  Sondrestrom  is  at 
a  higher  geomagnetic  latitude  than  the  other  two;  Sondrestrom  is  under 
the  oval  near  noon  and  in  the  polar  cap  from  dusk  until  dawn. 

In  addition  to  this  unique  property  during  MITHRAS  for  LT-UT 
studies,  the  observing  period  was  characterized  by  its  proximity  to 
solar-cycle  maximum  and  its  overlap  with  the  brief  lifetime  of  the 
Dynamics  Explorer-2  (DE-2)  satellite.  The  opportunity  to  obtain  simul¬ 
taneous  and  coordinated  data  from  three  radars  widely  spaced  around  the 
auroral  zone  was  the  driving  reason  for  the  timing  of  the  MITHRAS  exper¬ 
iment.  The  uniqueness  of  the  other  two  characteristics  only  became 
apparent  after  the  fact.  The  possibility  of  overlapping  with  DE-2  was 


considered  fortunate,  and  its  observations  were  coordinated  with  those 
of  the  radars.  However,  we  anticipated  that  the  satellite  would  remain 
in  orbit  longer  than  it  did.  Being  near  the  solar-cycle  maximum  was 
very  fortunate,  as  well.  Because  of  the  high  electron  densities  during 
this  period,  the  radar  data  set  is  of  very  high  quality.  Thus,  these 
data  were  unique  in  more  ways  than  were  originally  expected. 


The  radar  measurements  were  supplemented  by  those  from  STARE  (the 
Scandinavian  Jwin  Auroral  Eadar  Experiment),  the  Fabry-Perot  interfero¬ 
meter  network,  the  global  magnetometer  network,  the  DE-1  and  -2  satel¬ 
lites,  the  NOAA-6  and  -7  satellites,  the  IMP-8,  ISEE-1  and  -3  space¬ 
craft,  and  the  mid-  and  low-latitude  incoherent-scatter  radars.  They 
were  further  extended  with  older  data  from  Chatanika  and,  recently,  with 
data  from  Sondrestrom.  These  data  were  necessary  to  determine  back¬ 
ground  or  average  conditions  and  to  extend  the  coverage  into  the  polar 
cap. 


The  second  phase  of  MITHRAS  concerns  data  reduction  and  analysis. 
This  document  is  a  report  of  work  involving  SRI  personnel,  alone  and  in 
collaboration  with  others,  over  a  three-year  period  on  this  second 
phase.  Part  of  the  work  consisted  of  collecting,  reducing,  and  display¬ 
ing  data,  as  well  as  the  development  of  new  analysis  procedures.  These 
activities  are  discussed  in  Section  II.  The  other  part  Is  the  research, 
which  included  the  analysis  and  interpretation  of  radar  and  other  data. 
It  also  includes  the  comparison  of  data  with  the  results  of  model  calcu¬ 
lations.  Because  of  the  many  instruments  and  the  model  calculations, 
this  research  has  often  involved  collaboration  with  other  scientists, 
both  nationally  and  internationally.  This  research  is  reviewed  in 
Section  III.  The  initiative  for  these  projects  has  either  come  from 
individual  researchers  or  has  been  the  outgrowth  of  three  workshops  held 
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to  discuss  the  project.  These  workshops  proved  very  valuable  for  inter-  ■ 


mental  and  theoretical  backgrounds  in  the  MITHRAS  data  set.  The  minutes 
of  these  meetings  were  presented  in  the  Interim  Progress  Report  fWlckwar 
et  al. .  1984a].  A  list  of  researchers  who  participated  in  the  workshops 
or  studies  are  listed  in  Appendix  B.  The  research  efforts  have  led  to  a 
considerable  number  of  presentations,  reports,  and  publications.  They 
are  listed  in  Appendix  C. 
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II  DATA  CONSIDERATIONS 


Some  unusual,  but  anticipated,  data-handling  problems  have  arisen 
from  the  effort  to  bring  together  the  large  body  of  data  from  many 
ground-based  and  spacecraft  sensors  required  for  Project  MITHRAS.  A 
highly  adaptable  tape  format  had  to  be  developed  to  exchange  data  among 
the  MITHRAS  participants.  To  explore  the  large  amount  of  data  for  good 
situations  for  further  study,  it  has  been  necessary  to  develop  an 
extensive  set  of  display  programs,  mostly  in  color.  These  have  been 
used  to  make  an  overview  of  the  data  set  and  have  been  used  in  the 
scientific  analyses.  To  the  extent  possible,  all  the  radar  and 
correlative  data  have  been  brought  together  in  a  data  library  and 
catalog  at  SRI  International.  To  take  full  advantage  of  these  MITHRAS 
data,  we  have  had  to  extend  the  analysis  capabilities  considerably  in 
the  F  region  where  most  of  the  data  were  acquired.  To  look  for  the 
regular  patterns  of  variation,  we  have  developed  binning  and  averaging 
programs  to  use  on  large  bodies  of  data. 

A.  Data  Exchange 

In  a  project  that  involves  many  sources  of  data,  one  of  the  most 
important  considerations  is  to  be  able  to  exchange  and  manipulate  the 
data  easily.  While  this  comment  appears  obvious,  the  task  of  imple¬ 
menting  it  is  far  from  trivial.  Since  the  beginning  of  the  project,  we 
have  worked  closely  with  the  URSI  Incoherent-Scatter  Working  Group  and 
the  Incoherent-Scatter  Data  Base,  which  has  been  set  up  by  the  National 
Science  Foundation  at  the  National  Center  for  Atmospheric  Research 
(NCAR)  in  Boulder,  Colorado.  Indeed,  the  MITHRAS  work  on  data  exchange 
became  the  beginning  of  the  formal  international  data  base;  the  MITHRAS 
data  exchange  format  [de  la  Beau.lardlere  et  al..  1982]  that  we  developed 
under  this  contract  was  the  model  for  the  NCAR  data  base.  The  format 


has  now  evolved  and  been  adopted  for  the  NCAR  data  base,  as  well  as  by 
the  Europeans  for  the  EISCAT  data. 


Because  of  MITHRAS,  we  have  played  a  very  active  role  in  data 
exchanges.  It  has  also  meant  that  we  have  been  able  to  obtain  and 
readily  plot  the  data  from  the  other  radars  and  from  a  number  of  other 
instruments  with  which  we  need  to  compare  measurements.  It  has  also 
facilitated  exchange  of  data  with  theoreticians  and  modelers.  We  expect 
that  in  the  future  this  format  will  be  adopted  by  an  even  wider  segment 
of  the  aeronomy  community. 

B.  Data  Display 

The  radars  produce  an  enormous  amount  of  information.  This  is 
apparent  when  we  realize  that  the  measurements  can  be  reduced  and 
analyzed  to  produce  many  physical  parameters  as  a  function  of  altitude, 
latitude,  and  time.  The  problem  of  displaying  the  data  is  further 
exacerbated  because  some  of  the  parameters,  such  as  the  electron 
density,  can  vary  in  magnitude  over  such  a  very  wide  range  that  they  are 
hard  to  display. 

Perhaps,  if  we  knew  precisely  for  what  we  were  looking,  this  would 
not  matter.  But  this  is  an  observational  science  and  we  are  still  very 
much  exploring  the  high-latitude  region.  We,  therefore,  needed  to 
examine  as  many  of  the  parameters  from  as  much  of  the  data  as  possible. 
Indeed,  It  is  in  doing  so  that  we  have  found  a  number  of  unexpected 
results — for  example,  the  great  variation  of  nighttime  F-region  electron 
density  at  the  different  longitudes,  very  high  electron  temperatures 
near  the  convection  reversal  in  the  afternoon  convection  cell,  and,  on 
occasion,  coincidence  between  energetic  particle  precipitation  and 
F-region  electron-temperature  enhancement. 

To  help  in  displaying  the  data,  we  purchased,  as  part  of  Project 
MITHRAS,  a  color  terminal  and  SRI  purchased  a  camera.  We  developed  an 
extensive  series  of  color  display  programs.  The  color  termin' 1  is  an 
Advanced  Electronic  Design  (AED)  model  767  which  has  a  1024  x  1024  x  8 


bit  memory,  and  a  viewing  window  of  768  x  575  pixels.  Eight  bits  per 
pixel  allow  the  simultaneous  use  of  256  colors  from  a  palette  of  almost 
17  million  different  colors. 

The  color  camera  is  a  Dunn  631  system.  It  receives  the  three-color 
video  signals  from  the  AED  terminal  and  is  also  connected  to  the  main 
computer  by  a  serial  link. 

We  developed  and  documented  general-purpose  software  to  drive  these 
two  instruments  fLeeer .  1983].  This  user's  guide  describes  the  oper¬ 
ation  and  special  commands  for  the  AED  terminal  together  with  illustra¬ 
tive  examples. 

To  display  the  data  in  color,  the  first  practical  problem  was  to 
develop  a  color  scale.  With  the  ability  to  use  256  colors,  it  is 
possible  to  make  what  appears  to  be  a  continuous  scale.  However,  it  is 
easier  to  examine  data  where  discrete  steps  exist.  Therefore,  a  cri¬ 
terion  in  developing  a  color  scale  was  to  have  a  large  enough  color 
change  at  each  step  to  distinguish  between  colors.  Another  criterion 
was  to  have  a  "logical"  change  in  color  across  the  scale  so  that  we 
would  not  have  to  refer  back  to  the  scale  to  determine  if  a  color  step 
represented  an  increase  or  decrease. 

With  much  experimentation,  two  scales  were  developed.  The  first 
has  25  steps;  the  second  has  20.  The  25-step  scale  is  used  for  most  of 
the  data.  The  20-step  scale  is  used  for  velocity  data,  which  can  be 
positive  or  negative.  For  both  scales  it  was  necessary  to  vary  both  the 
hue  and  intensity  to  separate  the  colors.  The  "logical"  progression  in 
hue  is  in  the  following  order:  black,  blue,  cyan,  green,  yellow,  red, 
and  magenta.  The  subset  from  blue  to  red  is  the  usual  spectrum  obtained 
by  dispersing  white  light.  Furthermore,  with  the  exception  of  black, 
the  pattern  is  a  cyclic  progression  of  primary  colors  from  blue  to 
green,  green  to  red,  and  half  way  from  red  to  blue.  To  avoid  confusion 
between  the  highest  and  lowest  values,  the  cycle  is  not  completed  back 
to  blue.  The  intensity  settings  for  the  three  primary  colors  were  given 


^  [1984a].  They  may  be  useful  to  other  groups  beginning 


to  use  color. 


The  region  probed  by  the  radar  can  often  be  represented  as  a  prism, 
as  pictured  in  Figure  Il-l.  The  axes  are  altitude,  latitude,  and  time. 
During  the  MITHRAS  observation  campaign  three  main  experiment  modes  were 


used,  MITHRAS  1,  2,  and  3 


1982,  1984b].  The 


MITHRAS  1  and  3  experiments  fit  this  representation  well.  In  the 
MITHRAS  2  Experiment,  time  resolution  is  increased  at  the  expense  of 
latitude  coverage  and  the  separation  between  the  two  planes  defining  the 
triangle  becomes  very  small. 

In  terms  of  this  three-dimensional  solid,  it  is  apparent  that  we 
can  put  planes  through  the  solid  in  several  ways  to  get  useful  two- 
dimensional  representations  of  the  data.  We  have  developed  the  software 
to  do  this.  The  capabilities  are  illustrated  with  examples  using 
MITHRAS-3  data  from  Chatanifca  observations  on  27  January  1982. 

The  first  cut  through  the  data  is  a  latitude  cut.  Figures  II-2(a) 
through  II-2 (d)  show  the  electron  density,  electron  temperature,  ion 
temperature,  and  line-of-sight  ion  velocity,  respectively.  The  data 
were  acquired  during  an  elevation  scan  talcing  15  min.  The  data  were 
accumulated  for  55  s  and  reduced  to  the  basic  parameters  shown.  De¬ 
pending  on  the  quality  of  the  data  and  what  one  is  looking  for,  the 
accumulation  or  integration  time  could  be  shorter  or  longer.  The  axes 
are  altitude  and  invariant  latitude,  where  each  degree  of  invariant 
latitude  is  approximately  100  km  in  the  magnetic  north/south  direction 
at  the  earth's  surface.  Because  the  earth's  surface  has  been  flattened 
in  this  display,  the  radar  line-of-sights  appear  curved.  This  cut  is 
particularly  good  for  showing  the  relation  between  E-  and  F-region 
phenomena  and  for  showing  the  relative  latitudinal  locations  of 
structures  and  boundaries.  These  figures  show  E-region  ionization, 
indicating  energetic  particle  precipitation  between  64  and  67  A.  The 
elevated  F-region  electron  temperatures  in  coincidence  with  enhanced 
electron  density  also  indicate  soft  particle  precipitation  between 
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65  and  67  A.  Poleward  of  67  A,  the  hard  and  soft  precipitation  drop 


dramatically  and  Joule  heating  of  the  ions  sets  in.  In  these  repres¬ 


entations  structures  appear  vertical  when  the  inputs  are  aligned  along 


the  magnetic  field.  If  a  structure  were  not  vertical,  it  would  be  a 


warning  that  the  phenomenon  was  varying  rapidly  in  time  compared  to  our 


scan  rate. 


The  second  cut  through  the  solid  is  a  plane  along  the  magnetic 


field  direction,  with  altitude  and  time  as  the  two  coordinates. 


Figures  II-3(a)  through  II-3(c)  show  the  electron  density,  electron 


temperature,  and  ion  temperature.  Note  the  relation  between  E-  and 


F-region  parameters.  The  F-region  density  enhancement  before  midnight 


is  not  accompanied  by  energetic  particle  precipitation  that  creates  an 


E  region,  whereas  the  F-region  density  enhancement  after  midnight  is 
accompanied  by  energetic  particle  precipitation.  Similarly,  note  the 
covariation  or  lack,  of  covariation  among  parameters.  For  Instance,  no 
electron  temperature  increase  accompanies  the  density  increase  before 
midnight,  whereas  one  accompanies  the  Increase  after  midnight.  (The 
implication  is  that  the  density  increase  before  midnight  was  created 
elsewhere  and  transported  to  where  It  was  observed.  In  contrast,  the 
density  increase  after  midnight  was  being  created  locally.)  Another 
feature  that  stands  out  is  the  rapid  way  in  which  the  various  parameters 


turn  on  and  off  such  as  shortly  after  local  midnight  as  the  radar  moves 


from  one  region  to  another  or  as  particle  precipitation  and  convection 
turn  on  and  off.  That  some  of  these  changes  are  caused  by  crossing 


boundaries  will  be  more  obvious  in  the  next  cut. 


The  third  cut  through  the  solid  is  a  horizontal  plane  at  a  speci¬ 
fied  altitude.  The  axes  are  latitude  and  time.  The  plot  may  be  rec¬ 
tangular  or  curved  to  form  a  clock  dial.  Examples  of  electron  density 
at  250,  350,  and  450  km  are  shown  in  Figures  II-4(a)  through  II-4(c). 


Examples  of  electron  and  ion  temperature  at  350  km  are  shown  in 


Figures  II-5(a)  and  II-5(b),  respectively.  This  display  is  particularly 


good  for  showing  regions  and  boundaries.  For  instance,  an  ionization 


trough  develops  shortly  after  1400  LT  and  extends  through  the  midnight 
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sector.  The  way  the  latitude  of  the  trough  boundary  varies  in  time  can 

be  compared  to  the  convection  velocities.  Figure  II-6,  and  the  solar 

zenith  angles,  Figure  II-7,  all  presented  in  exactly  the  same  format. 

It  is  readily  apparent  that  the  trough  is  a  magnetospheric  related 

phenomenon,  i.e.,  a  phenomenon  related  to  convection,  not  to  solar 

insolation.  Indeed,  the  trough  appears  in  the  region  of  sunward 

convection  in  the  evening  convection  cell.  Poleward  of  the  trough, 

there  is  considerable  patchy  ionization.  Between  1900  and  2000  LT,  the 

elevated  electron  temperatures  strongly  suggest  that  the  patches  arise 

from  precipitation.  After  that,  the  low  temperatures  suggest  that  the 

patches  were  produced  elsewhere  and  convected  into  the  field  of  view. 

They  could  have  been  produced  just  poleward  of  the  observations  or  all 

the  way  across  the  polar  cap  near  noon.  After  midnight,  the  radar  moves 

under  the  region  of  sunward  convection  in  the  morning  cell.  It 

encounters  both  Joule  heating  (elevated  T  )  and  particle  precipitation 

(coincidence  in  the  occurrence  of  elevated  N  and  T  ). 

e  e 

We  have  presented  examples  of  the  three  major  ways  of  surveying  the 
data,  which  correspond  to  the  software  packages  that  were  developed 
under  MITHRAS.  In  the  examples,  we  have  shown  the  four  basic  parameters 
derived  from  the  radar  measurements.  Most  of  the  derived  parameters  can 
be  similarly  displayed.  In  addition,  by  using  the  exchange  format 
discussed  above,  we  can  obtain  and  display  data  from  other  radars  in  the 
same  way.  As  an  example,  In  Figures  II-8(a)  through  II-8(c),  we  show 
electron  densities  from  Chatanika,  Millstone  Hill,  and  EISCAT  for 
18  November  1981  fde  la  Beaulardlere  et  al..  1985a].  Major  differences 
in  the  diurnal  pattern  are  readily  apparent.  They  are  discussed  in  the 
next  section  in  terms  of  transport  across  the  polar  cap.  It  is  also 
apparent  that  the  high  nighttime  densities  in  the  F  region  at  Chatanika 
and  EISCAT  are  very  structured. 

In  addition  to  these  presentations,  we  have  added  the  ability  to 
include  some  correlative  data,  such  as  IMF,  on  the  same  plot.  An 
example  of  convection  velocities  and  IMF  is  presented  in  Section  III. 

We  have  also  added  the  ability  to  bin  and  plot  in  the  usual  ways  a  large 
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ION-CONVECTION  VELOCITIES  AT  CHATANIKA  ON  27  JANUARY  1982. 
The  format  is  the  same  as  in  Figure  11-5.  The  conventions  for  the  vectors 
are  that  the  tail  is  at  the  time  of  the  measurement,  that  the  head  is  in  the 
direction  of  the  ion  velocity,  and  that  the  vector  is  only  plotted  if  the 
uncertainty  in  its  magnitude  is  less  than  400  m/s.  For  example,  at  0600  UT 
the  ion  drift  is  greater  than  1  km/s  toward  magnetic  west  at  67° A,  while  at 
1245  UT  the  ion  drift  is  approximately  2  km/s  toward  magnetic  east  at  67°A. 


sample  of  data  so  then  we  can  examine  average  behavior  instead  of  parti¬ 
cular  cases.  An  example  of  the  variation  of  convection  velocity  with 
IMF  By  direction  is  presented  in  de  la  Beaulardlere  et  al.  [1985b; 

1986]. 

In  conclusion,  software  was  devised  to  display  the  data  in  three 
phases  from  the  four-dimensional  volume  formed  by  the  space  and  time 
coordinates.  This  software  has  provided  us  with  a  powerful  tool  for 
examining  the  various  high-latitude  phenomena. 


C.  Data  Analysis 

Before  MITHRAS,  our  research  efforts  concentrated  on  the  relatively 
easy-to-analyze  E-region  data  from  the  Chatanika  radar.  They  were 
principally  used  to  learn  about  energetic  electron  precipitation  and 
such  electrodynamic  parameters  as  conductivities  and  currents.  Although 
Chatanika  provided  much  F-region  data  as  well,  the  only  parameter  rou¬ 
tinely  used  was  the  easy-to-derive  convection  velocity.  However,  the 
fact  that  there  are  only  small  amounts  of  E-region  data  from  the  EISCAT 
experiments,  and  that  there  are  none  from  Millstone  Hill  in  the  auroral 
region  has  prompted  us  to  develop  analysis  procedures  for  using  the 
other  F-region  parameters  for  studying  auroral  related  phenomena:  soft- 
particle  precipitation,  downward  heat  flux,  electron-energy-loss  rates, 
optical  emissions,  and  meridional  neutral  wind. 

When  low-energy  electrons,  those  having  energies  less  than  400  eV, 
precipitate  into  the  atmosphere,  they  penetrate  to  F-region  altitudes 
only.  In  doing  so,  they  ionize  the  neutral  gas  and  heat  the  ambient 
electrons.  Because  of  the  long  chemical  lifetimes,  on  the  order  of  one 
to  several  hours,  F-reglon  ionization  can  be  transported  great  distances 
(Section  III.B).  Therefore,  increases  in  density  need  not  indicate 
particle  precipitation.  However,  the  electron  gas  heats  up  and  cools  in 
the  order  of  a  few  seconds  to  a  few  tens  of  seconds,  which  suggests  that 
electron  temperature  is  a  better  parameter  for  indicating  soft-particle 
precipitation. 
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However,  the  electron  temperature  alone  is  not  a  good  indicator  of 


energy  deposition  by  either  soft  particles  or  other  sources.  If  the 
electron  density  is  low,  a  given  energy  input  may  increase  the  electron 
temperature  greatly.  If  the  density  is  high,  that  same  energy  input 
will  raise  the  temperature  far  less.  Thus,  when  we  examine  electron 
densities  and  temperatures,  we  find  that,  if  two  regions  (at  the  same 
altitude)  have  the  same  electron  density,  and  one  has  a  higher  temper¬ 
ature,  then  more  energy  is  deposited  into  that  region.  Similarly,  if 
one  region  has  both  a  higher  density  and  a  higher  temperature,  then 
greater  energy  is  deposited  into  that  region.  However,  if  one  region 
has  a  lower  density  and  a  higher  temperature  than  another,  it  is  im¬ 
possible  to  say  from  a  qualitative  argument  alone  into  which  region  the 
greater  energy  is  deposited. 

Another  difficulty  with  a  qualitative  analysis  is  that  it  does  not 
tell  which  of  several  possible  mechanisms,  solar  EUV,  particle  precip¬ 
itation,  or  heat  conduction,  contributes  to  the  energy  deposition.  Most 
of  these  problems  can  be  solved  by  the  proper  quantitative  analysis  of 
the  data.  In  conjunction  with  the  Sondrestrom  radar  project,  we  have 
developed  the  software  to  make  the  above  analysis  more  quantitative. 

We  have  developed  least-squares  fitting  procedures  to  smooth  the 
electron-  and  ion-temperature  profiles  rKofman  and  Wichwar.  1984],  We 
have  also  devised  new  methods  for  estimating  the  neutral  temperature 
(and  its  asymptotic  value,  the  exospheric  temperature)  rAlcavde  et  al.. 
1983a;  Wlckwar  and  Kofman.  1984].  Using  these  profiles  and  a  model  of 
the  neutral  atmosphere,  the  energy  loss  rates  from  electrons  to  ions  and 
neutrals  can  be  calculated.  As  part  of  that  calculation,  the  downward 
heat  flux  into  the  electron  gas  near  500-fcm  altitude  can  also  be  derived 
fRofman  and  Wlchwar.  1984].  In  addition,  the  emission  at  6300  A  from 
atomic  oxygen  and  at  5200  A  from  atomic  nitrogen  can  be  calculated  for 
thermal  excitation  and  dissociative  recombination. 

Thus  we  can  determine  the  total  energy  loss  rate,  which  has  to  be 
equal  to  the  net  energy  input  to  the  electrons  from  all  sources,  local 
and  otherwise.  This  is  very  helpful  for  determining  the  location  of 


soft-particle  precipitation.  Although  most  of  the  energy  of  precip¬ 
itating  particles  may  not  be  carried  by  soft  particles,  they  are  im¬ 
portant  for  creating  F-region  ionization,  existing  F-region  optical 
emissions,  producing  metastable  species  that  can  be  transported  great 
distances,  and  for  affecting  the  thermospheric  wind  system.  We  can  also 
determine  the  energy  delivered  to  the  electron  gas  from  above  by  heat 
conduction.  Although  it  is  not  clear  whether  the  energy  is  deposited 
Just  above  our  maximum  altitude  by  very  soft  particles  or  whether  it  is 
conducted  from  far  out  in  the  magnetosphere,  this  is  a  new  capability 
that  we  have  never  had  before  in  the  auroral  region. 

We  have  significantly  improved  in  the  procedure  with  which  we  can 
derive  the  horizontal  neutral  wind  in  the  magnetic  north-south  dir¬ 
ection.  Our  ability  to  calculate  this  parameter  as  of  the  beginning  of 
MITHRAS  was  discussed  in  Wlcfcwar  et  al.  [1984c].  Since  then,  we  have 
improved  the  actual  calculation,  added  the  ability  to  use  the  Jacchia 
1971  model  and  to  vary  the  assumed  exospheric  temperature,  and  developed 
procedures  with  which  to  find  the  derivatives  of  the  temperatures,  which 
are  needed  to  determine  the  diffusion  velocity  properly.  Our  merid¬ 
ional-neutral-wind  calculations  have  now  been  developed  to  the  point 
where  they  are  a  very  useful  tool  for  ionospheric  and  thermospheric 
studies. 

The  derivation  of  analytic  temperature  profiles  are  useful  for 
additional  F-reglon  analyses.  We  can  now  calculate  the  excitation  of 
the  0(1D)  state  of  atomic  oxygen,  which  leads  to  emission  at  630  nm, 
from  dissociate  recombination  and  from  collisions  with  thermal  electrons 
[Wlckwar  and  Kofman.  1984]. 

We  have,  in  addition,  been  able  to  improve  the  derivation  of  the 
meridional  neutral  wind  f Wlckwar .  1984].  Instead  of  estimates,  the 
calculation  now  uses  the  electron-  and  ion-temperature  gradients  in  the 
diffusion  calculations. 


In  conclusion,  we  have  thus  considerably  extended  and  improved  our 
ability  to  use  the  F-region  data  from  the  radar  observations  to  learn 
about  the  physical  processes  occurring  in  the  F  region. 
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III  RESEARCH 


Appendix  C  lists  a  large  number  of  publications  and  presentations 
involving  MITHRAS.  These  are  the  result  of  new  observations,  efforts  to 
develop. and  validate  new  analysis  procedures,  studies  of  coincident 
radar  and  satellite  observations,  and  steps  in  an  ordered  progression 
from  simple  to  complex  comparisons  with  theory.  In  this  section,  we 
place  these  works  in  the  larger  context  that  we  outlined  in  Section  I, 
thereby  giving  a  "logical"  review  of  the  research. 

We  start  with  the  work,  that  pertains  to  the  magnetosphere,  namely, 
the  convection  pattern  and  the  precipitation.  From  there,  we  turn  to 
the  ionosphere,  i.e.,  the  F-region  electron  density  and  the  factors 
affecting  it,  the  electron  temperature,  and  the  ion  temperature.  Fin¬ 
ally,  we  discuss  the  neutral  atmosphere,  in  particular,  the  neutral  wind 
and  temperature. 

A.  Magnetosphere 

Two  major  manifestations  of  the  coupling  of  the  magnetosphere  to 
the  ionosphere  and  neutral  atmosphere  can  be  detected  with  the  radars. 
The  first  is  the  electric  field  that  is  mapped  down  the  field  lines  from 
the  magnetosphere  to  the  upper  atmosphere;  the  second  is  the  precipita¬ 
tion  of  electrons  and  ions  from  the  magnetosphere. 

i.  convection 

As  stated  above,  convection  plays  a  key  role  in  the  coupling  of  the 
magnetosphere  to  the  ionosphere  and  neutral  atmosphere.  Above  about 
180  km  both  ions  and  electrons  convect  in  the  £  x  2  direction.  A  number 
of  effects  occur  as  a  result  of  collisions  between  ions  and  neutrals. 
Momentum  is  transferred  to  the  neutrals,  dragging  them  along  with  the 
ions.  In  the  E  region,  the  ion  motion  is  impeded  by  collisions  with 
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neutrals.  Therefore,  the  ions  and  electrons  do  not  travel  at  the  same 
velocity  and  auroral  electrojet  currents  are  generated.  In  the  E  and 
F  regions,  the  ion-neutral  collisions  produce  Joule  heating.  The  energy 
from  Joule  heating  can  be  very  large.  We  will  look,  at  some  of  the 
effects  in  a  later  subsection. 

Our  interest  here  is  on  convection  itself,  not  its  consequences. 

We  will  describe  a  number  of  MITHRAS  studies  that  were  designed  to 
determine  the  convection  pattern,  what  affects  it,  and  how  it  evolves. 
Some  of  these  relied  on  simultaneous,  but  widely  separated,  observations 
by  the  high-latitude  radars;  others  depended  on  observations  from  the 
highest  possible  latitudes,  from  Sondrestrom. 

a.  Global -Convection  Pattern 

The  plasma  convection  pattern  in  the  polar  regions  is,  to  a  first 
approximation,  in  the  form  of  two  convection  cells;  the  flow  is  anti- 
sunward  across  the  polar  cap  and  sunward  along  the  auroral  oval  [Fig¬ 
ures  III-l(a)  to  (d)].  However,  considerable  uncertainty  and  contro¬ 
versy  still  exist  concerning  Important  aspects  of  plasma  convection. 

What  is  the  exact  shape  of  the  convection  pattern?  What  is  the  relative 
importance  of  the  dawn  and  dusk  cells,  what  are  the  precise  shapes  of 
the  cells,  and  is  the  overall  pattern  rotated  toward  early  hours?  Is 
this  pattern  substantially  modified  during  substorms  or  is  it  simply 
intensified?  Recently,  several  researchers  have  suggested  that  the 
pattern  may  be  more  complex,  consisting  of  three,  four,  or  even  five 
cells  intricately  distributed  over  the  high-latitude  region  rzanettl  et 
alt,  1984;  Potemra  et  al..  1984;  Relff  and  Burch.  1985;  Burgh  et  al^, 
1985].  These  more  complex  patterns  are  reproduced  in  Figure  III-l(e) 
through  (n)  for  several  IMF  orientations,  along  with  the  simpler  two 
cell  patterns  of  Frlls-Chrlstensen  et  al.  [1985]. 


It  is  well  established  that  the  IMF  strongly  influences  the  global 
convection  pattern.  However,  in  spite  of  a  large  number  of  obser¬ 
vations — ground  based  as  well  as  spaceborne — the  precise  role  of  each  of 
the  three  components  of  the  IMF  is  still  not  well  known.  The  role  of 


the  By  component  has  recently  received  considerable  attention  both  from 
a  theoretical  standpoint  f Crooker .  1979;  Banks  et  al . .  1984;  Relff  and 
Burch .  1985;  Lyons .  1985],  and  an  experimental  standpoint  [ Heells .  1984; 
Potemra  et  al..  1979;  Potemra  et  al..  1984;  Zanettl  et  al. .  1984;  Burch 
et  al. .  1985;  Jorgensen  et  al . .  1984;  Clauer  et  al . .  1984;  McCormac  and 
Smith.  1984;  Frlls-Chrlstensen  et  al..  1985].  However,  the  uncertainty 
about  how  the  convection  is  affected  by  IMF  B  Is  still  large.  This 
point  is  illustrated  by  the  variety  of  patterns  that  have  been  proposed 
in  three  recent  publications  on  the  subject,  which  are  used  as  sources 
for  Figures  III-l(e)  through  (n). 

One  reason  for  such  uncertainty  is  that  it  is  not  possible  to 
measure  the  plasma  convection  everywhere  at  the  same  time.  Satellite 
measurements  cover  a  wide  range  of  latitude,  but  only  along  one  line. 
Ground-based  measurements  with  the  incoherent-scatter  radars  have  the 
advantage  of  providing  time  variations  of  the  convection,  but  the  field 
of  view  is  small  compared  to  the  overall  convection  pattern.  At 
present,  the  most  effective  way  to  determine  the  global  pattern  is  to 
combine  data  from  several  radars  with  data  from  satellites.  Such  an 
analysis  was  done  in  collaboration  with  Heelis  f Heells  et  al. .  1983], 
Vector  velocity  measurements  from  DE-2  were  combined  with  MITHRAS  radar 
data.  Instantaneous  "snapshots"  of  the  polar  convection  were  derived  at 
the  time  of  each  northern  polar  pass  during  a  period  of  fairly  stable 
conditions  (25  October  1981).  One  of  these  snapshots  is  reproduced  in 
Figure  III-2,  which  shows  the  electrostatic  potential  distribution.  In 
the  F-region,  the  plasma  flows  along  these  equipotentials .  We  see  from 
Figure  III-2  that  the  afternoon  convection  cell  was  roughly  circular  and 
of  much  larger  extent  than  the  morning,  crescent-shaped  cell.  During 
these  observations  the  IMF  B z  was  negative  and  B  positive.  This 
orientation  corresponds  to  the  conditions  for  the  patterns  of  Figures 
III-l(a),  (c),  (e),  and  (g)  and  the  observed  convection  matches  the 
common  portions  of  patterns  (c)  and  (e)  fairly  well. 
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Further  studies  of  the  instantaneous  convection  pattern  are  needed. 


particularly  for  other  IMF  conditions  and  for  a  denser  set  of  observa¬ 
tions  such  as  is  obtained  using  Millstone  Hill,  Sondrestrom,  and  EISCAT. 


b.  Effect  of  the  IMF  By  Component 


Although  little  doubt  exists  that  the  IMF  By  component  strongly 


influences  the  whole  convection  pattern,  most  studies  have  concentrated 
on  the  daytime  pattern.  In  fact,  all  the  patterns  shown  in  Figure  III-l 
are  based  mostly  on  daytime  observations.  In  contrast,  little  attention 


has  been  given  to  the  By  effect  on  the  nightside  of  the  convection  since 


the  early  observations  of  Heotmer  [1972]  and  Mozer  et  al.  [1974]. 


In  three  separate  studies,  the  Sondrestrom  radar  24-hour  runs  have 


been  analyzed  to  determine  how  By  affects  the  nighttime  convection  [ de 


la  Beaulardlere  et  al..  1985b, c;  1986].  In  the  first  study,  four  events 
on  23  and  24  April  1983  were  examined.  In  the  other  studies,  two  years 
of  observations  from  1983  and  1984  were  examined,  both  individually  and 
statistically.  Velocity  vectors  from  the  first  study  are  shown  in 
Figure  III-3.  These  days  were  selected  because  they  are  representative 


of  features  observed  on  many  days.  The  IMF  B  and  B  vectors  (GSM 

z  y 


coordinates)  are  displayed  in  the  inside  circle.  Red  is  for  westward 


velocities  and  positive  B^;  blue  is  for  eastward  velocities  and  negative 


By.  B^  is  positive  toward  the  center  of  the  circle.  Averaged  vector 
velocities  are  shown  in  Figure  III-4.  These  figures  show  that  the  B 
component  of  the  IMF  strongly  influences  the  nightside  polar  convection. 
Its  effect  is  quite  complex.  The  convection  for  one  orientation  of  B 
is  not  the  mirror  image  of  the  other  orientation.  This  asymmetry  is  in 
disagreement  with  the  schematic  patterns  proposed  by  Potemra  et  al. 
[1979;  1984]  and  Relff  and  Burch  [1985],  but  it  is  present  in  the 
empirical  models  by  Heppner  [1984]  and  Fr lls-Chrlstensen  et  al.  [1985] 


(Figure  III-l).  A  positive  By  seems  to  organize  the  velocities  such 


that,  at  all  local  times,  they  are  predominantly  westward  within  the 
radar  field  of  view.  Figure  III-3(b)  offers  a  clear  example  of  this; 
the  radar  appears  to  remain  under  a  large  westward  vortex  for  24  hours. 
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Between  dusk,  and  midnight,  on  23  and  24  April  1983,  sunward  flow  is 
observed  within  the  polar  cap  [Figure  III-3(a)].  In  the  midnight  and 
dawn  sectors,  when  By  is  negative,  the  plasma  velocities  often  appear 
random  [18  January,  Figure  III-3(d)].  For  large  negative  By,  the  after¬ 
noon  cell  appears  shifted  toward  early  hours  such  that  large  southward 
velocities  are  observed  about  three  hours  before  midnight  [23  July, 
Figure  III-3(c)].  The  only  times  when  the  predominant  velocity  com¬ 
ponent  is  southward  are  in  the  premidnight  hours,  when  By  is  large  and 
negative. 

Thus,  the  B  effect  of  the  IMF  on  convection  is  very  clear  and  very 
dramatic.  This  avenue  of  MITHRAS  research  should  be  continued  and 
extended.  The  data  exist  to  do  so.  The  lower  latitude  data  from  Chat- 
anika,  pertaining  mostly  to  sunward  convection  and  the  region  around  the 
Harang  discontinuity,  should  be  examined. 

The  effects  of  other  factors  should  also  be  determined.  The  first, 
and  most  obvious,  is  season.  The  seasonal  dependence  of  currents  (iono¬ 
spheric  as  well  as  field-aligned  currents)  is  well  established,  but 
little  work  has  been  done  on  seasonal  dependence  of  polar  cap  convection 
velocities  or,  equivalently,  electric  fields.  Other  factors  Include  the 
various  activity  Indices.  Such  parameters  would  be  either  magnetic 
activity  indices  such  as  Kp,  AE,  AU,  or  global-energy  indices  such  as 
the  auroral  oval  radius  [Gussenhoven  et  al..  1983]  or  the  auroral  preci- 


the  auroral  oval  radius 
pitation  energy  index  [I 


1982].  The  various  energy-related 


parameters  of  the  solar  wind,  such  as  B  ,  VB  ,  or  e  ,  should  also  be 

z  z 

correlated  with  the  observed  electric  fields. 


Another  method  for  determining  the  large-scale  convection  pattern 
is  used  by  R.  Schunk,  J.  Sojka,  and  C.  Rasmussen  at  Utah  State  Uni¬ 
versity  (USU).  The  pattern  is  derived  by  adjusting  the  parameters  of  an 
analytic  expression  for  a  two-cell  representation,  such  as  those  of 
Volland  [1978]  or  Heel is  [1984],  until  a  good  fit  is  achieved.  In  a 
collaborative  study  between  USU  and  SRI,  Rasmussen  et  al.  [1986]  have 


v, 
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applied  the  technique  to  the  MITHRAS  data  of  27  and  28  June  1981  from 
Chatanlka,  Millstone  Hill,  and  STARE.  Despite  small  disturbances  (or 
substorms),  they  were  able  to  fit  one  pattern  remarkably  well  to  the 
whole  set  of  velocity  data.  This  pattern  Is  reproduced  In  Figure  III-5. 
It  Is  in  a  frame  of  reference  that  does  not  rotate  with  the  earth.  The 
validity  of  this  pattern  is  strongly  supported  by  the  fact  that  the 
derived  F-reglon  electron  densities  calculated  with  the  USU  global 


ionospheric  model,  which  depend  greatly  on  the  convection  pattern,  agree 


closely  with  the  observations.  This  Ng  comparison  will  be  discussed 


later. 


This  indirect  technique  for  inferring  the  large-scale  convection 
pattern  is  very  effective  and  more  such  studies  should  be  pursued.  We 
started  a  study  for  17  and  18  November  1981.  The  constraints  are 
greater  for  this  period  because  we  have  vector  ion  velocities  from  all 
three  MITHRAS  radars.  In  addition,  we  have  the  meridional  component  of 
the  thermospheric  neutral  wind  from  both  Chatanlka  and  EISCAT,  as  well 
as  two  lower-latitude  radars — St.  Santin  and  Arecibo.  This  period  is  a 
particularly  valuable  one  to  model  because  it  has  already  been  examined 
in  detail  in  two  separate  works  [Caudal  et  al..  1984;  de  la  Beaulardlere 
<?t  all  .  1985a]. 


Determination  of  the  Electrostatic  Potential 
Distribution  from  Magnetometers 


The  magnetic  disturbances  measured  from  a  worldwide  network  of 
magnetometers  can,  in  principle,  be  used  to  deduce  the  global  convection 
fKamlde  et  al..  1981;  Kamlde  and  Richmond.  1982].  Several  inversion 
methods  have  been  used  to  deduce  the  electrostatic  potential  from  the 
magnetic  perturbation  vectors  fFrlls-Chrlstensen  et  al..  1985  and  refer¬ 
ences  therein].  If  it  could  be  established  that  this  method  can  reli¬ 
ably  determine  the  global  convection  pattern,  it  would  be  very  powerful 
because  the  electric  potential  could  be  deduced  globally  and  with  good 
temporal  resolution.  However,  this  technique  is  severely  limited  by  the 
fact  that  the  results  strongly  depend  on  the  ionospheric  conductivities. 
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FIGURE  IH-5  SELECTED  PLASMA  STREAM  LINES  PREDICTED  FOR  27  AND  28  JUNE  1981  AS  A 
RESULT  OF  FITTING  CONVECTION  VELOCITIES  AT  CHATANIKA.  MILLSTONE  HILL, 
AND  STARE.  The  direction  of  motion  is  indicated  by  arrows;  one-hour  intervals  are 
indicated  by  dots.  The  coordinate  system  is  magnetic  local  time  versus  magnetic  latitude; 
the  frame  of  reference  is  quasi-inertial. 


Model  conductivities  based  on  statistical  averages  have  to  be  used 
because  of  the  paucity  of  measurements.  Despite  the  great  number  of 
magnetometers  distributed  over  the  globe,  problems  can  also  arise 

because  the  data  have  to  be  interpolated  between  the  magnetometer  | 

stations.  This  interpolation  can  introduce  significant  errors.  There-  \ 

fore,  the  question  still  remains  as  to  how  meaningful  the  results  of  i 

this  technique  are. 

The  multiradar  experiments  can  provide  very  useful  observations  for  ' 

improving  various  aspects  of  the  inversion  technique.  Work  was  started  j 

on  this  collaborative  project  using  MITHRAS  data  in  two  ways:  first  to 
provide  a  check,  on  the  inferred  plasma  convection,  and  second,  to  pro¬ 
vide  measurements  of  conductivities  that  can  be  incorporated  into  the 
calculations.  Particle  measurements  from  the  NOAA,  DMSP,  P78-1,  DE-2, 
and  HILAT  satellites  were  used  to  provide  estimates  of  the  conduct¬ 
ivities.  The  days  selected  were  25  October  1981,  26  October  1981, 

11  November  1981,  8  December  1981,  and  18  to  19  January  1984.  So  far, 

Y.  Kamide,  S.-I.  Akasofu,  and  B.-H.  Ahn  (Kyoto  Sangyo  University,  Uni¬ 
versity  of  Alaska,  and  Kyangpook  National  University,  respectively)  have 
collected,  calibrated,  digitized,  rotated,  and  verified  the  magnetometer 
data  from  over  100  sites  for  most  of  these  days.  A  plot  of  stacked 
magnetic  disturbances  from  about  60  stations  is  shown  in  Figure  III-6 
for  11  November  1981.  A  first  attempt  at  the  derivation  of  the  poten-  j 

tial  function  was  done  before  the  end  of  this  contract.  As  a  result,  i 

A.  D.  Richmond  (NCAR)  has  modified  the  inversion  technique  so  as  to 
include  the  measured  conductivities  and  electric  fields  whenever  and 
wherever  they  exist. 

e.  Electric  Fields  at  Subauroral  Latitude 

We  have  studied  subauroral  electric  fields  that  have  been  observed 
in  the  afternoon  or  evening  sector  during  several  MITHRAS  experiments. 

In  particular,  on  18  November  1981  f Caudal  et  al. .  1984]  and  11  November 
1981  ( Senior.  1982;  Senior  et  al. .  1986].  These  electric  fields  give 


rise  to  large  sunward  ion  drifts  that  occur  equatorward  of  the  auroral 


precipitation.  These  fields  occur  because  in  the  afternoon  sector,  the 
equatorward  boundary  for  proton  precipitation  is  not  collocated  with 
that  for  electron  precipitation  but  is  equatorward  of  it.  This  effect 
has  been  explained  in  terms  of  the  crosstail  electric  field  in  the 
magnetosphere  causing  a  sunward  displacement  of  the  inner  edge  of  the 
ring  current  rschield  et  al. .  1969;  Vasvllunas.  1972].  The  Region-2 


..  1969; 


Birkeland  currents  flow  between  these  two  boundaries.  They  close  in  the 
ionosphere,  in  a  region  of  reduced  conductivity.  As  a  consequence  of 
the  low  conductivity,  large  subauroral  electric  fields  are  established 
r Smiddv  et  al. .  1977].  Verification  of  this  explanation  requires  a 
careful  study  of  boundaries  for  proton  and  electron  precipitation  in 


relation  to  the  convection  pattern. 


Analysis  of  the  observations  from  11  November 


.,  1986] 


was  mostly  performed  while  Senior  was  at  SRI,  working  on  the  MITHRAS 
project.  This  particular  case  was  selected  because  there  were  coinci¬ 
dent  DE-2  and  NOAA-6  satellite  passes.  The  satellite  measurements 
provide  parameters  that  are  essential  to  complement  the  Chatanika  data: 
field-aligned  currents  as  well  as  proton  and  electron  precipitation 
spectra.  Figure  III-7,  from  Senior  et  al.  [1986],  illustrates  these 
observations.  The  figure  depicts  the  equatorward  boundary  of  the  dif¬ 
fuse  aurora  and  of  the  electric  field,  measured  by  the  Chatanika  radar. 
The  Region-1  and  Region-2  Birkeland  currents,  measured  by  DE-2  are  also 
shown. 
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Substorms  are  the  manifestations  of  energy  transfer  from  the  mag¬ 
netosphere  to  the  ionosphere.  This  transfer  can  be  an  explosive  release 
of  energy  that  has  been  previously  accumulated  in  the  tail  f McPherron . 
1979;  Rostoker  et  al..  1980]  or,  alternatively,  can  be  directly  driven 
from  the  solar  wind  energy  fAkosofu.  1980].  But  the  relative  importance 
of  these  two  components  is  not  known.  Whether  it  is  a  driven  system  or 
a  storing-unloading  of  energy,  this  energy  transfer  is  in  fact  a  complex 
superposition  of  a  variety  of  processes.  A  summary  of  a  substorm 
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FIGURE  III-7  DETAILED  CHATANIKA  AND  DE-2  OBSERVATIONS  IN  THE 
EVENING  SECTOR  AURORAL  ZONE  ON  11  NOVEMBER  1981. 

The  location  and  alignment  of  the  boundary  blob  with  the  convection 
and  diffuse  aurora  boundaries  strongly  suggest  that  it  is  related  to  the 
equatorward  edge  of  the  central  plasma  sheet. 


sequence  is  shown  in  Figure  III-8  I  NASA .  1984].  It  follows  from  this 
figure  that  a  complete  description  of  substorms  has  to  rely  on  a  variety 
of  observations  in  the  solar  wind,  magnetosphere,  and  ionosphere.  The 
difficulty  in  malting  the  appropriate  measurements  coincidently,  and  at 
enough  points  within  these  regions,  has  hampered  the  understanding  of 
how  these  processes  operate.  The  MITHRAS  observations  provide  a  very 
important  set  of  data  for  these  studies. 

The  advantage  of  the  radar  data  is  that  the  radars  provide  contin¬ 
uous  coverage  at  several  locations  and  that  they  can  provide  the  para¬ 
meters  that  are  needed  to  define  both  the  electrodynamic  and  precipi¬ 
tation  aspect  of  the  energy  transfer.  For  example,  radars  can  determine 
the  electric  fields  and  the  ionospheric  and  Birkeland  currents,  as  well 
as  the  energy  dissipation  parameters  such  as  Joule  heating  and  particle 
precipitation. 

Several  MITHRAS  studies  examined  the  convection  pattern  behavior  at 
two  or  three  radar  sites.  In  the  first  of  these  rde  la  Beaulardlere  et 
al, .  1983],  the  convection  velocities  at  Chatanika,  Millstone  Hill  and 
STARE  were  examined  during  five  periods  of  geomagnetic  disturbances  on 
27  to  28  June  1981.  At  the  same  local  time,  independent  of  longitude, 
similar  behavior  occurred.  In  the  dawn  and  duslt  sectors,  the  convection 
Intensified;  in  the  midnight  and  noon  sectors  it  appeared  to  reverse. 
When  Chatanika  and  EISCAT  were  in  the  dawn  and  dusk,  sectors,  respect¬ 
ively,  the  convection  intensified  at  both  locations,  almost  simultan¬ 
eously,  at  the  onset  of  a  magnetic  disturbance  or  substorm. 

The  analysis  of  a  series  of  magnetic  disturbances  was  also  carried 
out  for  the  experiment  on  27  January  1982  r Fontaine  et  al. .  1983;  de  la 
Beaulardlere  et  al..  1984a].  Again,  an  intensification  of  the  convec¬ 
tion  pattern  was  found  in  the  dusk,  and  dawn  sectors.  In  the  noon  per¬ 
iod,  the  convection  region  was  poleward  of  the  region  probed  by  the 
Chatanika  radar.  However,  significant  precipitation  of  energetic  part¬ 
icles  was  observed  in  the  noon  sector  above  Chatanika  at  the  time  of  the 
magnetic  disturbance.  Precipitation  by  low-energy  particles,  which 
deposit  their  energy  in  the  F  region,  was  examined  on  this  day  and  found 
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FIGURE  III-8  A  SUMMARY  OF  THE  SUBSTORM  SEQUENCE  SHOWING  THE  ROLES  OF  THE 
MAGNETOTAIL  AND  THE  IONOSPHERE  AND  ALSO  SHOWING  THE  ASPECTS 
OF  SUBSTORMS  THAT  ARE  REGARDED  AS  DRIVEN  AND  UNLOADING 
PROCESSES 


to  occur  in  conjunction  with  the  magnetic  disturbances  in  the  dusk,  and 
dawn  sectors.  This  was  the  first  attempt  at  using  a  new  method  to  infer 
the  existence  of  low-energy  precipitation. 


Caudal  et  al.  [1984]  performed  a  careful  comparison  of  convection 
velocities  from  three  MITHRAS  radars  for  18  November  1981.  The  con¬ 
vection  was  quite  steady  during  this  experiment.  On  this  stable  back¬ 
ground.  perturbations  were  observed.  When  EISCAT  was  under  the  after¬ 
noon  cell  and  Chatanika  under  the  morning  cell,  the  velocity  variations 
were  correlated. 

These  studies  show  that  during  substorms  sunward  convection  is 
increased  in  the  dawn  and  dusk  regions  and  often  velocity  reversals  are 
observed  near  noon  and  midnight.  It  is  too  early  to  say  whether  this 
behavior  is  due  to  an  increase  in  convection  velocities,  an  expansion  of 
the  convection  cells,  or  a  reconfiguration  of  the  pattern. 


Later  studies  related  the  convection  changes  to  changes  in  the  IMF. 
It  was  seen  that,  within  about  10  min  following  a  southward  turning  of 
the  IMF,  the  electric  field  intensifies  in  the  dawn  and  dusk  sectors  [de_ 
la  Beaulardiere  et  al..  1983,  1985b, c;  Clauer  et  al..  1984].  An  example 
of  timing  of  events  associated  with  a  change  in  the  IMF  orientation  is 
shown  in  Figure  III-9,  reproduced  from  de  la  Beaulardiere  et  al. 

[1985c].  Fifteen  minutes  after  the  IMF  turned  southward,  the  convection 
intensified  and,  to  the  north  of  Sondrestrom,  changed  direction.  Then, 
after  another  15  min,  a  sharp  onset  of  a  negative  bay  was  observed  in 
the  midnight  sector.  The  onset  of  a  negative  bay  at  the  radar  local 
time  occurred  much  later,  one  hour  after  the  IMF  orientation  changed. 
These  and  other  radar  observations  indicate  that  the  convection  inten¬ 
sifies  before  the  onset  of  a  substorm  expansive  phase.  But,  is  this 
intensification  part  of  the  substorm  growth  phase?  Is  it  simply  an 
intensification  of  an  existing  pattern  or  does  it  correspond  to  a  slg- 
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FIGURE  III-9  SEQUENCE  OF  EVENTS  FOLLOWING  IMF  CHANGE 
AT  2100  UT  ON  23  APRIL  1983 
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Thus,  a  number  of  case  studies  have  shown  the  potential  for  using 
the  radars  to  learn  about  substorm  processes  and  have  provided  infor¬ 


mation  about  convection  velocities  and  particle  precipitation.  The 
radars  are  a  very  powerful  tool  for  substorm  studies.  They  enable 
convection  and  particle  precipitation  to  be  examined  simultaneously  over 
a  range  of  latitudes.  Because  of  these  capabilities,  this  line  of 
research  can  continue  to  provide  significant  advances. 

2.  Precipitation 

A  second  major  aspect  of  coupling  between  the  magnetosphere,  on  the 
one  hand,  and  the  ionosphere  and  neutral  atmosphere,  on  the  other,  is 
the  precipitation  of  energetic  particles.  This  precipitation  leads  to 
auroral  emissions;  produces  D-,  E-,  and  F-region  ionization;  can  modify 
the  electric  fields;  and  deposits  a  substantial  amount  of  energy  in  the 
upper  atmosphere.  Indeed,  the  energy  deposition  can  be  very  significant 
in  local  areas,  compared  to  that  which  the  overhead  sun  would  deposit. 
Thus  it  can  have  large  effects  in  the  upper  atmosphere. 

Our  interest  in  this  section  is  on  the  precipitation  itself,  not 
its  consequences.  The  interpretation  of  the  MITHRAS  data  has  revealed 
where  various  types  of  precipitation  occur,  under  what  conditions,  and 
at  what  energies. 

The  radars  have  proved  well  suited  for  this  study.  Chatanifca, 
Sondrestrom,  and  EISCAT  are  sensitive  to  the  precipitation  of  electrons 
with  energy  between  a  few  tens  of  eV  and  100  fceV,  because  the  electron 
density  profile  can  be  measured  from  the  F  region  down  to  the  D  region. 
When  combined  with  photometers  or  spacecraft,  they  provide  information 
on  proton  precipitation.  In  addition,  they  enable  the  precipitation  to 
be  examined  in  the  larger  context  of  the  convection  pattern.  This  is 
particularly  important  because  precipitation  and  convection  are  closely 
linked,  although  in  the  past  they  have  almost  always  been  studied 
separately. 
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One  approach  to  examining  particle  precipitation  is  to  relate  where 


one  type  of  precipitation  is  located  with  respect  to  another.  This 
approach  was  followed  by  Wlnnlneham  et  al.  [1975;  1978]  with  satellite 
data.  They  were  able  to  relate  their  observations  to  the  known  particle 
populations  in  the  magnetosphere.  They  identified  two  principal  types 
of  precipitation:  boundary-  and  central-plasma-sheet  precipitation  (BPS 
and  CPS,  respectively).  An  extension  of  this  approach  is  to  relate  the 
particle  precipitation  to  the  convection  pattern.  An  analysis  of  this 
sort  was  performed  for  the  cleft  region  by  Relff  et  al.  [1978]. 

The  relative  locations  of  various  types  of  precipitation  and  of 
convection  can  be  deduced  from  radar  observations.  Initial  efforts  were 
made  with  Chatanika  data  using  ion  velocities  and  E-region  densities  in 
the  evening  sector  by  Robinson  et  al.  [1982]  and  in  the  early  morning 
sector  by  Senior  et  al.  [1982]  and  Kamlde  et  al.  [1984].  The  first 
paper  showed  the  diffuse  aurora  from  the  CPS  in  the  sunward  convection 
region  and  auroral  arcs  from  the  BPS  near  the  convection  reversal.  The 
second  paper  showed  that  arcs  were  present  in  the  sunward  convection 
region.  However,  most  of  the  time  these  data  were  limited  to  the 
regions  of  sunward  convection  on  the  night  side  of  the  earth. 

For  MITHRAS,  this  research  has  been  extended  by  including  data  from 
Sondrestrom.  Because  this  radar  is  farther  north,  the  convection  rever¬ 
sal  can  be  seen  more  often.  Under  this  contract,  the  early  afternoon 
sector  has  been  examined  by  Wlckwar  and  Kofman  [1984],  and  Kofman  and 
Wlckwar  [1984]  using  new  procedures  to  determine  the  presence  of  low- 
energy  electrons.  Considerable  soft-  and  hard-particle  precipitation 
has  been  found  in  the  vicinity  of  the  convection  reversal  in  the  after¬ 
noon  convection  cells.  The  precipitation  occurs  in  the  region  where  the 
gradient  of  the  current  would  be  negative  and  it  occurs  almost  entirely 
in  the  region  of  sunward  convection.  These  particles  probably  originate 
in  the  BPS.  At  the  BPS,  and  equatorward  of  it,  there  is  often  a  diffuse 
auroral  E  layer,  which  probably  arises  from  precipitation  from  the  CPS. 


b.  Very  Energetic  Particles 

At  the  other  end  of  the  energy  spectrum,  the  MITHRAS  data  contain 
information  on  more  energetic  particles,  those  that  penetrate  below 
95  km — often  to  75  lun.  Figure  III-10  shows  the  resultant  D-region 
densities  observed  on  a  quiet  and  on  an  active  day.  Eighteen  days  have 
been  examined  in  order  to  study  the  behavior  of  energetic  electron 
precipitation  and  of  D-region  electron  density.  Under  quiet  conditions, 
the  densities  vary  slowly  in  time  with  a  low-altitude  cutoff  between 
85  and  90  lun  near  noon.  In  winter,  the  altitude  cutoff  appears  similar 
near  noon,  but  the  densities  are  lower.  Often  energetic  auroral  pre¬ 
cipitation  occurs  near  midnight  or  near  dusk  or  dawn.  Near  dusk,  there 
is  a  tendency  for  a  distinct  layer  to  appear  near  95  km.  At  dusk  and 
dawn,  the  radar  should  have  been  equatorward  of  the  diffuse  aurora  and 
significant  convection;  hence,  inside  the  plasmasphere.  Significant 
energetic  precipitation  was  observed  on  some  days  near  noon.  Thus,  it 
appears  that  energetic  precipitation  originates  from  three  separate 
source  regions. 

On  some  days,  the  energetic  precipitation  occurs  for  limited 
periods,  but  on  others  it  is  present  for  many  hours.  When  this  ener¬ 
getic  precipitation  occurs,  it  appears  uniform  over  the  extent  of  the 
region  probed  by  the  radar  (50  to  75  km). 

B.  Ionosphere 

In  this  subsection,  we  discuss  studies  involving  the  effects  of 
convection  and  particle  precipitation  on  the  electron  density,  the 
electron  temperature,  and  the  ion  temperature.  We  also  discuss  some  of 
the  physical  processes  that  affect  these  quantities. 

1.  Electron  Density 

The  F-reglon  electron  density  at  high  latitudes  varies  greatly  as  a 
result  of  several  competing  factors.  Ionization  is  produced  by  solar 
EUV  and  by  soft-particle  precipitation.  Plasma  may  be  convected  over 
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large  distances.  The  F  layer  can  shift  up  or  down  in  altitude  depending 
on  the  direction  of  the  meridional  component  of  the  neutral  wind,  the 
zonal  component  of  the  electric  field,  or  an  impressed  electric  field 
parallel  to  the  magnetic  field.  The  altitude  profiles  are  modified  when 
the  electrons  or  ions  are  significantly  heated  by  particle  precipitation 
or  Joule  heating,  respectively.  The  ion  chemistry  can  be  greatly  modi¬ 
fied  as  a  result  of  large  differential  velocities  between  ions  and 
neutrals,  which  hasten  the  change  from  atomic  to  molecular  ions,  or  as  a 
result  of  changes  in  the  neutral  composition  that  are  produced  by  Joule 
heating. 


At  high  latitudes,  we  tend  to  think,  that  most  of  the  variations  in 
F-reglon  electron  density  occur  at  night  as  a  result  of  particle  precip¬ 
itation.  Although  dramatic  changes  do  occur  then,  very  significant  ones 
also  occur  during  the  day,  which  are  more  difficult  to  explain.  Such 
decreases  have  been  seen  in  ionosonde  and  satellite  data  fProlss.  1980]. 
The  MITHRAS  data  contain  several  examples  when  the  daytime  density  is 
reduced  from  one  day  to  the  next  by  a  factor  of  four.  Figure  III-ll 
shows  a  very  significant  depletion  at  two  altitudes  during  the  daytime 
on  25  October  1981.  The  densities  are  at  their  usual  levels  the  day 
before  and  the  day  after. 

These  variations  are  large  enough  and  the  data  comprehensive  enough 
to  be  amenable  to  detailed  study.  The  hypothesis  is  that  some  processes 
linked  to  magnetic  activity  have  increased  the  concentration  ratio  of  N2 
to  0  at  F-region  altitudes,  with  the  result  that  the  daytime  electron- 
production  rate  is  reduced  and  the  loss  rate  increased.  This  would  have 
happened  if  there  had  been  significant  Joule  heating  poleward  of  the 
radar,  followed  by  equatorward  transport  of  the  neutrals.  The  chemical 
part  of  this  mechanism  is  essentially  the  same  one  used  to  account  for 
larger  electron  densities  in  winter  than  in  summer  (the  winter  anomaly) 
and  for  greater  electron  densities  at  solar  cycle  maximum  than  at  solar 
cycle  minimum.  The  phenomenon  and  the  mechanism  are  most  likely  related 
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to  the  much  smaller  density  increases  and  decreases  at  midlatitudes 
after  a  magnetic  storm. 

b.  Longitudinal  or  UT  Density  Variations 

Because  of  the  separation  of  the  geographic  and  geomagnetic  poles, 
the  effects  of  solar  EUV  and  magnetically  controlled  phenomena  occur  in 
regions  that  are  partially  displaced  from  one  another.  This  displace¬ 
ment  is  even  more  complex  because  the  auroral  oval  is  shifted  about 
3°  antlsunward  of  the  magnetic  pole  rGussenhoven  et  al..  1983].  These 
offsets  lead  to  a  complicated  relative  motion  of  the  auroral  zone  with 
respect  to  the  geographic  pole.  In  turn,  this  motion  leads  to  longi¬ 
tudinal  or  UT  variations  in  the  electron  density  behavior.  Such  vari¬ 
ations  have  been  predicted  rKnudsen.  1974;  So llta  et  al. .  1982]  and 
observed  rsato  and  Rourke.  1964]. 

One  of  the  important  MITHRAS  observations  was  that  in  the  middle  of 
the  night  in  fall  and  winter,  the  F-region  density  at  EISCAT  was  greater 
than  at  Chatanika  and  the  density  at  Chatanika  was  greater  than  at 
Millstone  Hill,  Figure  III-12.  These  longitudinal  variations  have  been 
examined  in  considerable  detail  by  de  la  Beaulardlere  et  al.  [1985a]. 

The  data  are  consistent  with  the  idea  that  ionization  produced  in  the 
vicinity  of  the  noon  auroral  oval  and  cleft  can  be  transported  across 
the  polar  cap  and  that  the  production  of  ionization  near  the  cleft  by 
solar  EUV  varies  because  of  the  offset  between  the  geomagnetic  and 
geographic  poles.  In  other  words,  the  relative  location  of  the  term¬ 
inator  and  auroral  oval  varies  with  longitude. 

Several  other  UT  variations  have  been  observed  in  the  MITHRAS  data 
set.  One  involves  the  rate  of  density  increase  at  dawn  and  decrease  at 
sunset.  At  Chatanika,  the  increase  is  slow  and  the  decrease  rapid.  At 
EISCAT,  the  reverse  is  true.  Another  is  the  noon  density,  which  is 
often  larger  at  Chatanika  than  at  EISCAT.  A  third  is  the  E-region 
electron  density,  with  the  result  that  the  conductivities  appear  to  be 
systematically  greater  at  EISCAT  than  at  Chatanika. 
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FIGURE  III- 1 2  MAGNETIC-LOCAL-TIME  DEPENDENCE  OF  THE  F-REGION  ELECTRON 

DENSITY  AT  MILLSTONE  HILL,  CHATANIKA,  AND  EISCAT  ON  18  AND 
19  NOVEMBER  1981  AND  15  AND  16  DECEMBER  1981.  The  densities  are 
from  an  altitude  of  500  km  and  an  invariant  latitude  of  approximately  64°. 
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There  are  several  large-scale  F-region  density  features  that  are 
closely  associated  with  the  convection  pattern.  As  such,  they  may  also 
have  a  longitudinal  dependence.  The  best  studied  is  the  trough 
[Muldrew.  1965],  an  electron  density  depletion  that  has  a  limited  lati¬ 
tudinal  extent.  There  can  be  multiple  depletions.  The  "main  trough"  or 
"midlatitude  trough"  is  in  the  evening  or  nighttime  sectors,  but  deple¬ 
tions  occur  at  any  local  time  [Grebowskv  et  al..  1983]. 

The  trough  has  been  examined  extensively  with  data  from  Millstone 
Hill  [Evans  et  al. ■  1983a,  1983b;  Holt  et  al. .  1983].  They  discuss 
three  major  factors  that  appear  to  cause  the  observed  depletion.  The 
first  is  the  stagnation  hypothesis  [Knudsen.  1974].  Because  sunward 
convecting  plasma  in  the  evening  sector  spends  a  long  time  in  darkness, 
the  normal  chemical  decay  has  a  long  time  to  operate.  The  second  is  the 
subauroral  electric  field  hypothesis  f Evans  et  al. .  1983a].  These 
fields,  as  already  discussed,  occur  in  the  evening  sector  and  lead  to  an 
increased  rate  for  the  charge  exchange  reaction  0+  +  N2  +  N0+  +  Ne 
f Banks  et  al..  1974;  Schunk  et  al..  1975],  which  leads  to  an  accelerated 
F-region  decay.  The  third  is  the  fossil  hypothesis  r Evans  et  al. . 
1983b].  It  is  similar  to  the  previous  hypothesis  in  that  it  depends  on 
increased  charge  exchange  followed  by  recombination  to  cause  the 
decrease  in  ionization.  It  differs  in  that  it  depends  on  large  electric 
fields  during  magnetic  disturbances  or  substorms  to  create  the 
depletion,  which  then  persists  and  corotates  once  the  disturbance  has 
subsided  and  the  region  of  increased  convection  has  moved  poleward. 

The  main  trough  is  on  the  equatorward  side  of  the  auroral  oval 
where  it  is  observed  in  the  region  of  sunward  convection  in  the  after¬ 
noon  and  evening  sectors.  But  the  effect  of  large  ion  velocities  on  the 
chemistry,  which  is  involved  in  the  second  and  third  hypotheses,  can 
occur  elsewhere.  Indeed,  large  reductions  in  density  are  often  seen  at 
Sondrestrom  poleward  of  the  auroral  oval,  in  regions  of  strong  anti- 
sunward  convection.  The  velocities  are  often  between  1  and  3  km/s. 
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We  have  carefully  examined  a  MITHRAS  data  set  in  which  the 
midlatitude  trough  extended  over  several  hours  of  local  time  (Fig¬ 
ure  III-13).  This  case  study  was  done  while  Senior  was  at  SRI.  In  it. 
Senior  et  al.  [1986]  found  that  the  trough  formation  was  consistent  with 
the  stagnation  model  developed  by  Spiro  et  al.  [1978]. 

In  the  same  study,  we  explained  the  formation  of  a  ridge  of  ioniz¬ 
ation,  which  was  observed  to  lay  parallel  to  the  main  trough,  along  its 
poleward  wall.  We  showed  that  this  ridge  of  ionization  (often  called 
the  boundary  blob)  could  not  arise  from  local  production,  but  was  the 
result  of  the  large-scale  convection. 

The  MITHRAS  observations  are  more  comprehensive  than  those  which 
were  previously  available.  The  data  from  Chatanika,  Sondrestrom,  or 
EISCAT  have  better  latitudinal  resolution  than  the  data  from  Millstone 
Hill,  which  were  used  for  the  early  radar  trough  studies.  The  data  from 
Chatanika  and  EISCAT  have  a  wide  local-time  separation.  Data  from 
Sondrestrom  are  from  a  higher  latitude,  enabling  the  trough  and  other 
structures  to  be  seen  at  earlier  times  as  well  as  enabling  the  deple¬ 
tions  formed  in  regions  of  antisunward  convection  to  be  explained.  More 
can  be  done  with  these  data  as  procedures  are  developed  for  a  detailed 
comparison  between  theory  and  observation.  The  status  of  these  compar¬ 
isons  is  discussed  in  the  next  section. 


The  last  aspect  of  our  study  of  ionospheric  densities  concerns 
collaborative  work  with  the  USU  group.  As  was  mentioned  in  an  earlier 
section,  we  just  performed  the  first  comparison  of  their  model  with  data 
from  several  radars  f Rasmussen  et  al. .  1986].  This  work  combines  the 


from  several  radars  I  Rasmussen  et  al. .  1986].  This  work  combines  the 
USU  model  ionosphere  with  data  from  27  June  1981,  principally  from  the 
Chatanika  and  Millstone  Hill  radars,  but  also  includes  data  from  the 
STARE  radar  and  the  NOAA-6  satellite.  Despite  the  occurrence  of  several 
magnetic  disturbances,  it  was  possible  to  use  a  time-independent  con¬ 
vection  pattern.  In  general,  good  agreement  was  obtained  between  obs¬ 
ervation  and  theory,  indicating  that  no  major  physical  or  chemical 
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lb)  ELECTRON  DENSITY  AT  350  km  IN  UNITS  OF  cm‘J 

HGURE  111-13  MAIN  ELECTRON  DENSITY  TROUGH  AT  CHATANIKA  ON 
11  NOVEMBER  1981  The  main  trough  moves  equatorward 
from  72  A  at  0130  UT  to  60  A  at  0700  UT  Several  other 
regions  of  density  depletion  occur  during  the  night 
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process  was  left  out  or  badly  calculated.  There  were,  however,  a  number 
of  smaller  discrepancies  that  should  be  looked  at  more  carefully  in 
future  studies. 

In  this  case  study,  the  convection  pattern  and  neutral  wind  pattern 
were  held  constant  for  24  hours.  Only  ion  velocities  and  electron 
densities  were  more  extensive  than  in  previous  calculations  in  that 
velocities  were  available  from  two  locations  all  the  time  and  from  three 
locations  some  of  the  time.  Similarly,  the  electron  densities  were 
compared  at  several  altitudes.  A  follow-up  study  was  started  for  the 
same  day  to  explore  electron  and  ion  temperature  behavior. 

From  the  experience  gained  in  these  studies,  the  next  step,  when 
funding  is  renewed,  will  be  to  allow  the  convection  pattern  to  change 
continuously  throughout  the  period  being  modeled.  This  will  allow  days 
with  more  variable  geomagnetic  activity  to  be  examined.  The  other 
change  that  will  have  to  be  tried  is  varying  the  neutral  wind  pattern  as 
a  function  of  time.  At  that  time,  it  will  also  become  necessary  to 
combine  the  ionospheric  model  with  the  thermospheric  global  circulation 
model  ( TGCM ) .  It  should  then  be  feasible  to  carefully  examine  such 
phenomena  as  the  daytime  depletion,  the  several  UT-longitudinal  vari¬ 
ations,  and  the  troughs.  The  combination  of  the  multiradar  measurements 
and  the  global  ionospheric  modeling  is  very  powerful,  but  the  process  is 
Just  beginning  and  the  ability  to  use  it  to  understand  more  complex 
situations  is  growing  incrementally. 

2-  Electron  Temperature 

Another  useful  parameter  to  study  in  order  to  gain  insight  into 
coupling  phenomena  is  the  electron  temperature.  This  parameter  has  the 
advantage  over  F-region  density  that  the  time  constants  for  heating  and 
cooling  are  on  the  order  of  a  few  seconds  instead  of  tens  of  minutes. 
Therefore,  it  is  a  good  indicator  of  energy  input  to  the  electrons. 
Generally,  this  energy  input  is  from  collisions  with  suprathermal  elec¬ 
trons.  These  electrons  may  be  photoelectrons  produced  by  solar  EUV  or 
UV,  or  secondary  electrons  produced  by  energetic  auroral  precipitation. 
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Another  source  of  energy  Input  is  heat  flux  down  the  magnetic  field. 

The  energy  source  may  arise  from  processes  far  out  in  the  magnetosphere 
or  from  auroral  particles  with  such  low  energy  that  they  penetrate  to 
Just  above  the  field  of  view  of  the  radar.  The  electrons  are  cooled  by 
collisions  with  ions  and  neutrals  as  well  as  by  heat  conduction.  Ulti¬ 
mately,  moat  of  the  energy  deposited  in  the  electron  gas  is  deposited  in 
the  neutral  gas. 


We  found  that  the  F-region  electron  temperatures  are  significantly 
o  o 

elevated  in  a  2  to  3  region  near  the  convection  reversal  in  the  after¬ 
noon  convection  cell.  We  have  seen  this  enhancement  near  noon  and  in 


the  early  afternoon  at  Sondrestrom 


1,  1984;  Wlckwar 


and  Kofman.  1984];  and  in  the  evening  and  near  midnight  at  Chatanika 
rFonthelm  et  al. .  1984;  Wlckwar  et  al..  1984d].  In  these  cases,  signi¬ 
ficant  energy  is  transferred  to  the  electrons.  This  energy  transfer  is 
readily  apparent  qualitatively  because  the  temperature  and  density  both 
increase.  It  has  been  examined  quantitatively  by  solving  the  electron- 
energy  equation  r Kofman  and  Wlckwar.  1984;  Wlckwar  et  al..  I984d]. 
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Under  steady-state  conditions,  the  rate  of  energy  input  to  the  electron 
gas  Qg  is  equal  to  the  rate  of  energy  loss  from  the  electron  gas  Lg  and 
the  rate  of  energy  loss  from  heat  conduction.  The  energy  input  can  be 
from  ionization  by  solar  EUV  and  UV  or  from  collisions  with  auroral 
primaries  and  secondaries.  The  energy  loss  can  be  from  collisions  by 
the  ambient  electrons  with  ions  and  neutrals.  The  net  energy  input  rate 
at  a  given  height  is  equal  to  L@.  The  radar  data  can  be  used  to  derive 
enough  parameters  to  calculate  L@  and  the  heat  conduction  term,  which  is 
closely  related  to  the  heat  flux. 

Under  extreme  conditions,  the  energy  deposited  between  180  and 

2 

500  km  has  reached  1  erg/cm  s.  If  it  were  from  soft  electrons,  they 
would  have  to  have  energies  below  400  eV  to  have  stopped  in  the 
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F  region.  Under  these  same  conditions,  the  downward  heat  flux  at  450-km 


altitude  was  0.3  erg/cm  s.  More  typical  values  were  an  order  of  magni¬ 


tude  smaller. 


In  MITHRAS  data  from  Chatanika  and  EISCAT,  we  have  found  that  there 


often  is  often  a  high  correlation  between  the  occurrence  of  energetic 


precipitation  producing  an  auroral  E  region  and  the  occurrence  of  ele¬ 


vated  electron  temperatures  in  the  F  region.  This  correlation  is 
apparent  in  the  data  from  the  two  parts  of  Figure  III-14.  Part  (a) 


shows  the  F-region  temperatures;  Part  (b)  shows  the  E-region  densities. 


We  have  started  to  interpret  these  observations  to  determine  under  what 


conditions  the  correlation  exists  and  to  determine  whether  the  elevated 


electron  temperatures  arise  because  of  soft  particle  precipitation  or 


heat  conduction. 


As  mentioned  above,  the  largest  electron  temperature  increases  have 
been  seen  near  the  convection  reversal  in  the  afternoon  cell.  Fig¬ 
ure  III-15  shows  the  extreme  examples.  The  first  is  from  near  noon  on 
24  April  1983  at  Sondrestrom  fWlckwar  and  Kofman.  1984].  The  second  is 


from  the  evening  on  5  March  1981  at  Chatanika 


.,  1984b]. 


In  the  second  example,  which  was  at  night,  an  extremely  bright  red 


aurora  covered  the  sky.  In  both  cases,  we  obtained  electron  temperature 
as  high  as  6000  K  at  500-km  altitude.  In  addition  to  the  high  temp¬ 
eratures,  we  demonstrated  that  large  downward  heat  fluxes  were  present 


on  both  days.  The  temperatures  and  fluxes  are  shown  in  Figure  III-15. 


We  are  now  trying  to  understand  why  such  large  temperature  in¬ 
creases  were  observed  on  these  two  days  and  not  on  other  days.  At  this 


point,  what  we  know  is  that  magnetic  activity  was  very  high  on  both 


days.  During  the  Sondrestrom  observations,  Kp  was  6-;  during  the 


Chatanika  observations,  a  large  magnetic  storm  was  in  progress. 


Other  questions  have  arisen  concerning  the  heat  flux.  Does  it 
arise  from  very  soft  particles  (about  10  to  20  eV)  stopping  at  altitudes 
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just  above  500  km,  or  does  it  arise  from  processes  deep  within  the 
magnetosphere?  Our  approach  to  these  questions  is  to  make  a  very 
detailed  comparison  of  the  Chatanika  data  and  one-dimensional  model 
calculations.  This  collaboration  has  been  ongoing  with  M.  Rees  and 
G.  Romick  (University  of  Alaska)  and  others. 

A  similar  question  is  being  addressed  on  several  MITHRAS  data  sets 
using  radar-determined  temperatures  in  the  topside  ionosphere  below 
550  km  and  DE-2  determined  temperatures  at  higher  altitudes.  The  work 
has  been  conducted  in  collaboration  principally  with  E.  Fonthelm  of 
Michigan  and  L.  Brace  of  NASA.  This  work  follows  a  study  using  data 
from  the  August  1972  magnetic  storm  ( Fonthelm  et  al. .  1978].  The 
present  work  is  aimed  at  learning  whether  an  anomalous  thermal  conduct¬ 
ivity  is  needed  to  explain  the  very  high  electron  temperatures  r Fonthelm 
et  al..  1984]. 


c.  Excitation  of  the  Red  Line  of  Atomic  Oxygen 

One  of  the  ways  that  hot  electrons  lose  their  energy  is  in  collis¬ 
ions  with  neutrals.  We  have  studied  the  possibility  that  hot  electrons 
could  excite  atomic  oxygen  to  the  state  followed  by  emission  at 
630.0  nm.  In  examining  the  radar  data  from  Sondrestrom,  we  found  that 
thermal  excitation  did  appear  capable  of  producing  several  kilo  Ray¬ 
leighs  of  red-line  emission  rwickwar  and  Kofman.  1984].  These  inten¬ 
sities  are  comparable  to  what  have  been  observed  in  cleft  auroras.  Our 
findings  are  important  because  we  found  that  the  electron  temperatures 
were  much  higher  than  assumed  by  previous  works  (Link  et  al. .  1983; 

Roble  and  Rees.  1977].  Indeed,  the  electrons  are  hot  enough  that  ther¬ 
mal  excitation  becomes  a  significant  excitation  mechanism. 

We  then  performed  the  same  calculations  with  data  from  5  March 
1981,  the  day  with  several  periods  of  intense  red  aurora.  We  again 
found  that  hot  electrons  could  account  lor  several  kR  of  red-line 
emission.  Temperature  and  heat  flux  prefiles  were  already  shown  in 
Figure  III-15.  Unlike  the  Sondrestrom  observations,  when  it  was  sunlit, 
the  Chatanlka  observations  were  at  night,  when  there  were  simultaneous 
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photometric  observations  by  researchers  from  the  University  of  Alaska 
and  the  Lockheed  Palo  Alto  Research  Laboratory.  The  optical  observa¬ 
tions  have  been  cross  calibrated  and  we  are  making  a  detailed  comparison 
of  Intensities  deduced  from  the  radar  data  and  those  that  were  meas- 


Preliminary  results  are  shown  in  Figure  III-16  for  data  obtained 
along  the  magnetic  field  line.  The  optical  data  (curve  labelled  "Photo¬ 
meter")  were  selected  from  very  rapid  meridian  scans  such  that  the 
fields  of  view  and  integration  times  were  as  closely  matched  as  possible 
to  the  radar  measurements.  This  ability  to  compare  radar  and  optical 
data  properly  is  an  important  analysis  capability.  We  calculated  the 
630.0-nm  intensities  from  two  mechanisms  (curve  labelled  "Radar"  in 
Figure  III-16 ) .  The  first  is  based  on  chemistry  rwickwar  et  al..  1974; 
Cogger  et  al..  1979];  the  second  is  based  on  collisions  with  hot  ambient 
electrons  fWickwar  and  Kofman.  1984].  When  auroral  activity  was  very 


electrons  IWickwar  and  Kofman.  1984].  When  auroral  activity  was  very 
low,  the  background  intensity  owing  to  the  recombination  of  ions 
accounts  for  much  of,  if  not  all,  the  observations.  The  largest  inten¬ 
sities  coincide  with  the  periods  of  greatest  electron  temperature,  when 
the  thermal  excitation  mechanism  accounts  for  the  observations.  Thus, 
we  have  confirmation  that  thermal  excitation  can  be  an  important  emis¬ 
sion  mechanism.  Under  intermediate  conditions,  i.e.,  with  lower  elec¬ 
tron  temperatures,  the  dominant  excitation  mechanism  is  presumably  from 
impact  of  the  atomic  oxygen  by  the  low-energy  precipitating  electrons 
and  their  secondaries. 

Because  of  the  importance  of  thermal  excitation,  the  ratio  of  the 
630.0-nm  intensity  to  that  of  other  auroral  emissions  is  not  as  good  a 
way  of  determining  the  precipitating  electron  energy  spectra  as  has  been 
assumed  for  the  last  decade.  This  finding  has  important  implications 
for  a  number  of  studies  because  the  ratio  of  intensities  at  6300  A  and 
4278  A  has  become  a  standard  technique  for  estimating  a  Maxwellian 
distribution  for  energetic  auroral  electrons. 
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3.  Ion  Temperatures 


The  ion  temperature  is  another  parameter  that  can  provide  insight 
into  coupling  phenomena.  As  was  the  case  for  the  electron  temperature, 
the  time  constants  for  heating  and  cooling  the  ions  are  the  order  of  a 
few  seconds;  therefore,  the  ion  temperature  is  a  good  indication  of 
energy  being  transferred  into  and  out  of  the  plasma.  Ions  are  heated 
above  the  neutral  temperature  by  Coulomb  collisions  with  hot  electrons. 
This  process  has  noticeable  effects  above  about  300  km  and  takes  place 
at  all  latitudes.  Another  source  of  heating  predominant  at  high  lati¬ 
tudes  is  frictional  or  Joule  heating.  It  arises  from  collisions  between 
ions  and  neutrals,  most  usually  as  the  ions  are  driven  through  the 
neutrals  by  convection  electric  fields.  The  heating  rate  is  propor¬ 
tional  to  the  square  of  the  difference  between  the  vector-ion  and 
vector-neutral  velocities.  This  process  has  noticeable  effects  above 
115  km.  Thus  enhanced  ion  temperatures  are  very  useful  for  Indicating 
when  Joule  heating  is  important. 

a.  Elevated  Ion  Temperatures  in  the  Morning  Sector 

We  analyzed  coincident  Millstone  Hill  and  Chatanika  observations  of 
ion  temperatures  and  electron  densities.  We  found  that,  consistently, 
the  highest  observed  ion  temperatures  were  in  the  early  morning  r Baron 
and  Wand.  1983].  This  result  is  illustrated  in  Figure  III-17.  It  shows 
that  for  a  given  convection  velocity,  the  ion  temperature  is  higher  in 
the  morning  than  in  the  evening  sectors. 

The  explanation  given  by  Baron  and  Wand  [1983]  is  that  in  the 
evening  sector,  the  convecting  ions  rapidly  drag  the  neutrals  along  with 
them  thereby  reducing  the  velocity  difference  and  the  Joule  heating.  In 
contrast,  in  the  morning  sector  the  velocity  difference  between  the  ions 
and  neutrals  persists  much  longer,  giving  rise  to  more  Joule  heating. 

We  attributed  the  difference  in  ion  drag  between  morning  and 
evening  to  larger  ion  densities  in  the  evening  sector.  This  study  was 
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[Source:  Baron  and  Wand.  1983] 

FIGURE  111-17  PREFERENTIAL  ION  HEATING  IN  THE  EARLY  MORNING  AT 
CHATANIKA  AND  MILLSTONE  HILL.  Although  the  range  of  ion 
velocities  is  similar  in  the  morning  and  evening,  the  ion  temperatures 
are  much  hotter  in  the  morning. 


pursued  by  Alcavde  et  al.  [1984].  Using  other  MITHRAS  data,  they  con¬ 
cluded  that  the  Coriolis  force  is  also  a  contributing  factor.  This 
early  morning  Joule  heating  probably  affects  the  neutral  atmosphere. 

The  energy  deposited  in  the  neutrals  could  raise  the  temperature,  and 
the  increased  temperatures  give  rise  to  pressure  changes  that  affect  the 
neutral  wind  circulation.  The  energy  deposition  can  also  lead  to  an 
upwelllng  that  changes  the  neutral  composition.  These  types  of  changes 
most  probably  produce  the  daytime  density  depletion  discussed  earlier. 

If  so,  they  should  also  lead  to  important  changes  equatorward  of  the 
auroral  region. 

b.  Comparison  with  Theory 

The  model  calculations  by  the  USU  researchers  [e.g.,  Rasmussen  et 
al. .  1986]  mentioned  in  earlier  sections  have  been  extended  to  include 
Joule  heating.  The  collaborative  work,  comparing  MITHRAS  data  and  model 
calculations,  was  started.  It  was  expected  that  the  results  would  help 
in  our  effort  to  analyze  the  complex  factors  that  affect  the  ion  and 
neutral  temperatures. 


c.  Neutral  Atmosphere 

As  mentioned  earlier,  the  magnetosphere  and  ionosphere  are  closely 
coupled  to  the  neutral  atmosphere.  Because  the  neutrals  dominate  at 
thermospheric  altitudes,  it  is  important  to  understand  the  behavior  of 
the  neutral  parameters  for  providing  background  conditions  and  for 
Identifying  the  largest  transfers  of  energy  and  momentum.  Convectlng 
ions  will  drag  the  neutrals  along,  significantly  changing  the  thermally 
driven  wind  pattern.  Particle  heating  and  Joule  heating  affect  the 
neutral  composition,  temperature,  and  winds.  These,  in  turn,  affect  the 
ion  and  neutral  chemistry.  In  the  following  subsections,  we  discuss 
studies  involving  the  neutral  wind  and  temperature. 


The  component  of  the  thermospheric  neutral  wind  in  the  magnetic 
meridian  is  derived  from  radar  measurements  parallel  to  the  magnetic 


field.  The  meridional  wind  can  be  derived  at  several  altitudes  in  the 
F  region,  throughout  the  day,  and  during  both  winter  and  summer.  The 
major  questions  concerning  the  thermospheric  wind  are  still  the  diurnal 
pattern,  how  it  is  affected  by  solar  variations  and  geomagnetic  activ¬ 
ity,  and  what  are  the  major  driving  forces. 

a.  Meridional  Wind  at  Chatanlka 

The  original  MITHRAS  study  of  the  meridional  wind  is  by  Wlckwar  et 
al.  [1984c].  Much  of  this  effort  went  into  developing  the  analysis 
procedure.  Basically,  the  ion  velocity  parallel  to  the  magnetic  field 
arises  from  neutrals  colliding  with  ions  and  from  ion-neutral  diffusion. 
To  derive  the  wind  from  the  observations  involves  calculating  the  ion- 
neutral  diffusion,  which  involves  the  derivatives  of  electron  density, 
electron  temperature,  and  ion  temperature  with  respect  to  altitude.  It 
also  Involves  a  careful  determination  of  the  neutral  atmosphere.  Good 
agreemt it  was  found  between  the  derived  wind  and  the  same  component  of 
the  wind  derived  from  optical  observations  with  a  Fabry-Perot  inter¬ 
ferometer. 

Chatanlka  data  from  seven  years  were  analyzed  and  compared  to 
theory.  The  average  equinox  pattern  was  established:  the  wind  was 
northward  during  the  day  at  about  65  m/s  and  southward  during  the  night 
at  about  175  m/s.  During  periods  of  increased  magnetic  activity,  the 
wind,  on  average,  became  more  southward  during  the  night  by  about 
100  m/s.  On  individual  days,  the  increase  was  as  much  as  400  m/s. 
Comparison  with  theory  showed  basic  agreement,  but  suggested  possible 
heating  poleward  of  Chatanika  during  the  morning. 


The  MITHRAS  observations  have  revealed  the  existence  of  anomalous 
behavior  in  the  nighttime  meridional  wind.  These  include  a  secondary 
equatorward  maximum  during  the  night  in  the  auroral  region  (Fig¬ 
ure  III-18)  and  an  abatement  at  high  latitudes  (Figure  III-19).  In  the 
top  portion  of  Figure  III-18,  the  wind  is  strongest  toward  the  south  at 
1200  UT,  or  0200  LT.  This  is  a  very  quiet  day,  as  indicated  by  a  north¬ 
ward  electric  field  that  is  always  less  than  10  mV/m.  The  wind  results 
are  typical  of  magnetically  quiet  days.  The  bottom  portion  of  the 
figure  shows  a  local  maximum  in  the  wind  between  1000  and  1200  UT, 
followed  by  a  much  stronger  maximum  near  1400  UT.  This  second  maximum 
coincides  in  time  with  the  strongest  southward  electric  field,  which  was 
50  mV/m.  This  second  maximum  is  typical  of  geomagnetlcally  active 
days. 

Figure  III-19  shows  the  meridional  wind  at  Sondrestrom  for  a  32-hr 
period  in  April  1983.  At  about  0200  UT  (2300  LT) ,  the  wind  is  toward 
the  north.  Usually  at  this  time  of  the  night  it  is  strongly  southward, 
but  on  some  days  there  is  an  abatement  and  on  this  night  a  reversal. 
Other  examples  of  this  behavior  have  been  presented  in  a  collaborative 
study  with  the  University  of  Michigan  f Meriwether  et  al..  1984]. 

We  believe  that  these  effects  are  indicative  of  the  influence  of 
magnetospheric  convection  on  the  neutrals.  These  effects  are  probably 
the  result  of  ion  drag  and  the  shape  of  the  convection  pattern  over  the 
polar  cap,  for  the  case  of  Sondrestrom  observations,  or  over  the  auroral 
zone  near  midnight,  for  the  case  of  Chatanifca  observations.  However, 
there  is  still  a  possibility  that  they  indicate  energy  deposition  such 
as  will  be  discussed  in  the  next  subsection. 

Our  early  research  emphasized  the  wind  determination  and  the  large- 
scale  meridional  wind  pattern  rwickwar  et  al..  1984c].  The  anomalous 
behavior,  just  discussed,  stands  out  in  comparison  with  this  pattern. 
Research  has  been  started  in  order  to  investigate  these  two  cases.  To 
do  so,  the  convection  pattern,  particle  precipitation  and  Joule  heating 
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FIGURE  III- 1 8 


SECOND  MAXIMUM  IN  THE  MERIDIONAL  WIND  NEAR  MID¬ 
NIGHT  UNDER  DISTURBED  CONDITIONS  AT  CHATANIKA. 
The  top  panel  shows  one  maximum  under  quiet  conditions.  The 
bottom  panel  shows  a  strong  second  enhancement  between  0200 
and  0530  LT  under  active  conditions. 
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have  been  observed  over  an  extended  latitude  region  on  either  side  of 
the  radar.  The  radar  data  have  been  supplemented  with  neutral  wind  data 
from  a  Fabry-Perot  interferometer  deep  in  the  polar  cap  at  Thule.  This 
last  effort  is  in  collaboration  with  the  University  of  Michigan. 


c.  The  Role  of  Particle  and  Joule  Heating 

Although  the  importance  of  ion  drag  rFuller-Rowell  et  al..  1981; 
Heppner  and  Miller.  1982;  Roble  et  al. .  1982;  Havs  et  al. .  1984;  Wlckwar 
et  al. .  1983]  in  determining  the  high-latitude  thermospheric  wind  has 
become  very  clear,  we  have  found  situations  at  Sondrestrom  that  strongly 
suggest  that  energy  input  from  both  particle  precipitation  and  Joule 
heating  affect  the  neutral  wind  rwicfcwar .  1984].  These  results  provide 
the  impetus  to  further  examine  the  daytime  and  nighttime  data.  As  in 
the  previous  subsection,  the  appropriate  observations  come  from  experi¬ 
ments  that  survey  a  wide  range  of  latitudes.  It  is  then  feasible  to 
correlate  the  behavior  of  the  overhead  wind  to  the  particle  and  Joule 
heating  poleward  and  equatorward  of  the  radar.  Sufficient  cases  exist 
so  that  a  combination  of  case  studies  and  a  statistical  analysis  would 
greatly  increase  our  Knowledge  about  these  effects. 


d.  B  Dependence  of  the  Neutral  Wind 

We  performed  a  statistical  analysis  on  all  the  data  from  a  two-year 
period  for  which  the  IMF  data  also  existed.  The  results  are  shown  in 
Figure  III-20.  This  figure  shows  the  average  neutral  wind  at 
Sondrestrom  for  two  orientations  of  the  IMF  By  component.  The 
differences  are  statistically  significant. 

As  noted  above,  the  wind  is  often  greatly  affected  by  ion  drag. 
Because  the  shape  of  the  convection  pattern  strongly  depends  on  the  IMF 
By  component,  it  is  not  very  surprising  that  the  wind  also  depends  on 
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FIGURE  III-20  By  DEPENDENCE  OF  THE  MERIDIONAL  WIND.  The  data  are  from 
Sondrestrom  from  the  period  March  1983  to  March  1985. 


We  performed  a  similar  statistical  analysis  to  look,  for  seasonal 
effects.  Because  most  of  the  Chatanika  data  were  from  the  equinoxes,  we 
again  relied  on  the  Sondrestrom  data  which  included  sufficient  obser¬ 
vations  from  the  solstices  to  make  this  analysis  possible.  The  results 
are  shown  in  Figure  III-21.  The  solid  curve  is  for  the  summer;  the 
dashed  curve  for  the  winter.  Between  the  early  morning  and  noon,  the 
wind  is  more  northward  in  the  summer  than  in  the  winter. 

This  difference  is  consistent  with  what  is  expected  from  a  stronger 
solar  input  in  summer  than  in  winter.  However,  the  fact  that  the  dif¬ 
ference  begins  so  early  in  the  morning  suggests  that  the  explanation  is 
more  complex.  Collaborative  work,  combining  observation  and  theory,  was 
started  in  order  to  seek  an  explanation  for  these  results. 

f .  Altitude  Gradients 

The  altitude  dependence  of  the  F-region  wind  had  not  been  studied 
before.  The  usual  working  assumption  was  that  the  wind  was  constant 
with  altitudes  above  about  200  km.  This  assumption  is  necessary  in 
order  to  Interpret  the  optical  data  because  the  emission  altitude  is  not 
well  known.  We  found  that  with  the  MITHRAS  data  this  assumption  is  true 
under  some  conditions,  but  not  others.  Figure  III-22  illustrates  this 
point.  It  shows  gradient  in  the  midnight  sector. 

Wlckwar  et  al.  [1984c]  had  previously  shown  that  gradients  tended 
to  appear  in  the  radar  data  during  the  nighttime  equatorward  maximum. 
However,  it  is  only  recently  that  the  analysis  program  has  been  improved 
to  the  point  that  we  are  sure  the  gradient  is  not  an  artifact  of  the 
calculations.  We  believe  that  the  gradient  is  closely  associated  with 
ion  drag  during  nighttime. 

In  going  from  200-  to  300-km  altitude,  the  electron  density  usually 
increases  and  the  neutral  density  always  decreases.  Therefore,  the 
ratio  of  ions  to  neutrals  increases  greatly  with  altitude,  with  the 
effect  that  ion  drag  is  more  effective  at  the  high  end  of  the  altitude 
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FIGURE  III-21  SEASONAL  DEPENDENCE  OF  THE  MERIDIONAL  WIND.  The  data  are  from 
Sondrestrom  from  the  winters  of  1983/84  and  1984/85  and  from  the  summers 
of  1983  and  1984.  Winter  and  summer  are  considered  to  be  three  months 
long,  centered  about  their  respective  solstices. 
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FIGURE  III-22  GRADIENT  IN  THE  MERIDIONAL  NEUTRAL  WIND  AT  CHATANIKA  ON 
18  NOVEMBER  1981.  The  curve  at  1002  to  1105  UT  (0002  to  0105  LT) 
shows  the  gradient;  the  other  curves  do  not.  The  error  bars  are  one  standard 
deviation  of  the  mean.  They  are  estimated  from  seven  3-minute  observations 
during  a  63-minute  interval. 


rang*.  The  nighttime  period  when  the  gradient  appears,  moreover,  coin¬ 
cides  with  when  the  peak  of  the  F  layer  is  at  its  highest  altitude  and 
when  Chatanika  is  under  the  portion  of  the  convection  pattern  with  the 
strongest  southward  ion  velocities. 

g.  Comparison  with  Theory 

In  addition  to  detailed  examinations  of  the  data  under  selected 
conditions,  another  approach  has  been  started  that  involves  an  in-depth 
analysis  of  one  or  two  days,  which  combine  as  many  data  as  possible  with 
extensive  modeling.  The  day  on  which  work  was  started  was  18  November 
1981.  Very  good  coverage  at  high  latitudes  exists  at  Chatanika  and 
EISCAT.  Good  coverage  at  middle  latitudes  also  exists  at  St.  Santin  and 
Arecibo.  At  high  latitudes.  Millstone  Hill  also  has  convection  and 
electron  density  Information. 

2.  Neutral  Temperature 

The  second  neutral  parameter  of  interest  is  the  neutral  temperature 
or  its  asymptotic  value,  the  exospheric  temperature.  This  parameter 
indicates  changes  in  the  neutral  atmosphere.  At  the  same  time,  a  sig¬ 
nificant  change  in  the  neutral  temperature  is  indicative  of  a  very  large 
energy  input  to  the  atmosphere.  It  can  lead  to  major  readjustments  in 
the  atmosphere  that  can  affect  the  chemistry  over  wide  areas  and  for 
extended  periods. 

a.  Horning  incxgftg.e 

As  mentioned  earlier,  we  found  that  large  ion  temperatures  often 
occur  in  the  early  morning.  They  are  attributed  to  Joule  heating.  This 
process  also  deposits  energy  in  the  neutrals.  The  question  is  whether 
the  neutrals  also  undergo  a  significant  temperature  increase.  We  have 
found  considerable  indirect  evidence  that  they  do.  Several  deter¬ 
minations  of  the  exospheric  temperature  from  the  radar  data  fWlckwar  et 
al. .  1984c;  Kelly  et  al. .  1983]  do  suggest  that  significant  increases 
occasionally  occur  in  the  morning  sector.  Figure  III-23  from  Wlckwar  et 
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FIGURE  III-23  EARLY  MORNING  ENHANCEMENT  OF  EXOSPHERIC 
TEMPERATURE  AT  CHATANIKA.  The  line  of  short 
dashes  in  both  panels  is  from  the  CIRA  [1972]  model 
atmosphere.  AST  is  Alaska  Standard  Time.  The  top 
panel  shows  temperature  variation  similar  to  that  of  CIRA. 
The  bottom  panel  shows  a  large  morning  enhancement 
relative  to  CIRA. 


al.  [1984c]  shows  three  Instances  of  this  behavior.  In  the  morning,  the 
exospheric  temperature  was  200  K  above  what  would  have  been  predicted  by 
an  empirical  neutral  model.  For  comparison,  we  show  that  the  temper¬ 
ature  follows  the  predicted  behavior  on  two  other  days. 

The  existence  of  these  early-morning  temperature  increases  are 
supported  by  other  observations.  A  few  optical  measurements  have  re¬ 
vealed  comparable  or  larger  increases  [Meriwether,  private  communi¬ 
cation,  1984].  Anomalous  increases  in  the  equatorward  meridional  wind 
in  the  early  morning,  which  have  already  been  discussed,  are  consistent 
with  significant  heating.  In  addition,  a  theoretical  effort  to  recon¬ 
cile  electric  fields  and  neutral  winds  at  Chatanika  had  to  include  such 
a  temperature  increase  in  the  early  morning  l~Havs  et  al. .  1979]. 

Thus,  it  is  likely  that  the  neutral  atmosphere  is  heated  sub¬ 
stantially  in  the  early  morning  and  that  the  mechanism  is  Joule  heating, 
but  the  case  has  not  yet  been  adequately  proved.  It  would  be  important 
to  verify  this  hypothesis. 

The  exospheric  temperature  is  derived  from  the  electron  and  ion 
temperatures  in  a  manner  similar  to  that  described  by  Bauer  et  al. 
[1970].  The  derivation  is  done  by  solving  the  ion  energy  equation  when 
it  is  possible  to  assume  that  the  only  energy  input  to  the  ions  is  from 
Coulomb  collisions  with  electrons.  However,  at  high  latitudes  the  ions 
are  heated  by  collisions  with  neutrals — Joule  heating.  Working  with 
groups  at  ONES  and  CNET,  we  have  attempted  to  correct  for  Joule  heating 
rAlcavde  et  al..  1983a].  However  the  correction  is  uncertain  because  of 
changes  in  ion  composition.  As  a  result,  the  safest  procedure  for 
determining  the  exospheric  temperature  is  to  consider  only  those  times 
and  places  where  the  Joule  heating  is  expected  to  be  negligible  rwickwar 
and  Kofman.  1984]. 

b.  Behavior  During  Decreasing  Phase  of  Solar  Cycle 

Another  study  involves  the  decrease  in  exospheric  temperature  from 
the  time  of  the  peak  of  the  solar  cycle  to  the  present,  i.e.,  starting 
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with  the  MITHSAS  data.  These  long-term  variations  are  Important  for  the 
behavior  of  the  neutral  atmosphere  and  Ionosphere.  We  have  found  that 
during  the  period  from  1981  to  1985,  the  exospheric  temperature  has 
dropped  by  about  700  K,  from  1400  K  to  700  K. 

During  this  same  period,  the  high-latltude  F  layer  has  had  roughly 
a  factor-of-four  decrease  in  maximum  electron  density  and  a  100-km 
decrease  in  the  altitude  of  the  layer.  These  changes  in  ionospheric 
behavior  are  closely,  though  not  exclusively,  tied  to  the  neutral 
atmosphere.  In  addition  to  these  general  reasons  for  determining  the 
exospheric  temperature  decrease,  preliminary  work  suggests  that  the 
changes  may  not  be  smooth.  For  example,  we  observed  a  sudden  decrease 
near  the  end  of  MITHEAS,  around  January  1982  rAlcavde  et  al..  1983b]. 
Whether  this  abrupt  decrease  was  a  perturbation  or  part  of  a  monotonic 
decrease  needs  to  be  examined  with  later  data. 


IV  DEVELOPMENT 


As  discussed  In  various  places  in  the  previous  sections,  knowledge 
of  the  solar  wind  and  interplanetary  magnetic  field  (IMF)  is  vital  for 
understanding  almost  any  auroral  zone  or  polar  cap  phenomena.  However, 
it  has  been  increasingly  difficult  to  measure  and  retrieve  the  solar 
wind  parameters. 

The  only  spacecraft  that  can  regularly  monitor  the  solar  wind  is 
currently  IMP-8.  In  the  past,  the  International  Sun  Earth  Explorer, 

I SEE- 3 ,  provided  continuous  measurements  of  solar  wind  parameters,  but 
in  June  1982,  its  orbit  was  changed  to  the  distant  magnetotail  and  then 
to  the  Giacovini-Zinner  comet.  The  ISEE-1  and  -2  and  the  AMPTE-IRM 
orbits  are  such  that  they  can  monitor  the  solar  wind  for  only  a  few 
weeks  each  year.  Furthermore,  ISEE-1  and  -2  are  expected  to  reenter 
during  1987.  The  WIND  spacecraft,  which  is  a  part  of  the  International 
Solar  Terrestrial  Program,  will  be  well  equipped  to  monitor  the  solar 
wind,  but  it  will  only  be  launched  in  the  nineties.  Therefore,  at  this 
time,  the  solar  wind  can  only  be  monitored  by  IMP-8. 

Unfortunately,  NASA  does  not  acquire  IMP-8  data  continuously 
because  its  telemetry  is  on  a  frequency  that  is  not  supported  by  the  new 
TDRS  network.  Recognizing  the  importance  of  the  solar  wind  data,  tele¬ 
metry  is  being  received  at  Wallops  Island  and,  very  recently,  at  the 
European  Space  Agency  Tracking  station  in  Belgium.  But,  to  extend  the 
coverage,  data  must  be  acquired  from  the  Pacific  sector.  For  this 
reason,  we  received  year-end  funds  to  modify  the  SRI  150-ft  antenna,  so 
that  it  could  acquire  the  IMP-8  telemetry,  thereby  extending  the  daily 
coverage  by  3  to  4  hours. 

The  antenna  was  built  by  SRI  in  1962  to  study  meteor-burst  communi¬ 
cations.  In  1964  and  1965,  a  computer-based  tracking  system  employing  a 


PDP-8  computer  was  added.  Other  improvements  were  also  added  for  study¬ 
ing  interplanetary  and  planetary  electron  densities  by  transmitting  to 
instruments  on  Pioneer  spacecraft  (Pioneer  5  through  9)  and  Mariner  5 
(Venus  mission). 

The  antenna  system  is  located  in  the  foothills  above  Stanford 
University  and  has  the  following  characteristics: 

•  Type — parabolic  reflector 

•  Diameter— 150  ft 

•  Design  frequency — VHF  through  UHF 

•  Mount — azimuth  elevation 

•  Gain— 32  dB  at  137  MHz 

•  Tracking — program  track  from  punched  paper  tape 

•  Location— 37. 402659°N.  -122.174139°E 

•  Altitude— 145.085  m. 

Work  was  performed  on  a  number  of  elements  in  the  antenna  system. 

In  particular,  a  new  feed  system  was  designed,  built,  and  installed. 
Repairs  were  made  on  the  old  PDP-8  computer  and  on  the  tracking  system. 
The  rotary  Joints  were  repaired.  The  antenna  system  was  tested.  The 
tests  showed  that  in  spite  of  the  intense  industrial  noise  from  Silicon 
Valley,  the  received  signal-to-noise  ratio  is  sufficient  to  acquire  the 
IMP-8  telemetry. 

Therefore,  it  has  been  determined  that  the  150-ft  antenna  and 


receiver  system  could  be  used  to  receive  the  IMP-8  telemetry.  To 
actually  use  the  system  reliably,  some  further  improvements  would  be 
needed,  as  would  funding  for  operations. 


V  CONCLUSIONS 


MITHRAS  has  proven  to  be  particularly  valuable:  not  only  did  it 
greatly  increase  our  scientific  understanding  of  the  upper  atmosphere  at 
high  latitudes,  but  it  also  increased  our  ability  as  a  group  of 
researchers  to  address  global  scientific  problems.  A  data  acquisition 
phase  was  well  thought  out  and  coordinated;  as  a  result,  a  valuable  and 
sizable  body  of  data  was  collected  under  a  variety  of  geophysical 
conditions.  Through  the  MITHRAS  workshops  and  presentations  at 
meetings,  many  researchers  participated  in  the  physical  interpretation 
of  the  observations.  In  addition,  arrangements  were  made  for  the  easy 
exchange  and  transfer  of  the  data.  As  a  result  of  these  efforts,  other 
data  sets  have  been  made  available  and  good  working  relationships  have 
been  established  among  several  experimenters  and  theoreticians.  Thus, 
MITHRAS  because  a  collaborative  effort  among  groups  who  work  with 
incoherent-scatter  radars,  as  well  as  with  a  large  number  of  other 
researchers.  Sufficient  support  has  been  provided  since  the  data 
acquisition  phase  to  allow  for  significant  research  activity  and  to 
enable  development  of  new  analysis  techniques. 

The  MITHRAS  project  had  an  important  role  in  the  establishment  of 
the  incoherent-scatter  data  base,  which  is  located  at  NCAR.  It  has  also 
served  as  a  model  for  a  number  of  smaller  multiradar  research  projects, 
such  as  GISMOS  (global  Incoherent-Scatter  Measurements  gf  Substorms)  and 
GTMS  (global  Thermospheric  Mapping  Study)* 

The  MITHRAS  experiment,  which  focused  on  an  extensive  set  of 
incoherent-scatter  data  from  the  high-latitude  region,  produced  a  number 
of  important  scientific  studies  that  increased  our  knowledge  about  the 
solar  wind,  magnetosphere,  ionosphere,  thermosphere,  and  the  inter¬ 
actions  among  these  regions: 
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•  The  best  view,  to  date,  has  been  obtained  of  the  instantaneous 
shape  of  the  convection  pattern  in  the  auroral  region  and  polar 
cap.  Instantaneous  snapshots  of  the  large-scale  convection  were 
obtained  by  combining  Millstone  Hill,  Chatanlka,  and  DE-2 
electric-field  measurements. 

•  A  series  of  case  studies  established  that  the  convection  inten¬ 
sifies  simultaneously  in  the  morning  and  afternoon  convection 
cells  during  a  substorm. 

•  A  case  study  of  the  convection  at  Chatanlka,  Millstone  Hill,  and 
El SCAT  showed  that  the  large-scale  pattern  could  remain  stable 
for  a  prolonged  period,  even  though  the  period  was  one  of  sus¬ 
tained  magnetic  activity. 

•  Average  patterns  of  plasma  convection  were  calculated  for  two 
different  orientations  of  the  IMF  By  component.  The  patterns 
obtained  showed  that  the  effect  of  the  By  component  is  stronger 
and  different  than  previously  thought. 

•  At  or  near  the  convection  reversal  in  the  afternoon  cell,  i.e., 
at  the  polar  cap  boundary,  there  can  be  large  energy  Inputs  to 

the  electron  gas  by  soft  particles  and  a  downward  heat  flux. 

•  The  winter  nighttime  F-region  density  strongly  depends  on  longi¬ 
tude;  this  dependency  was  attributed  to  cross-polar  cap  trans¬ 
port  of  ionization  and  to  the  fact  that  because  of  the  offset  of 
the  geographic  and  geomagnetic  poles,  the  solar  illumination  of 
the  large-scale  convection  pattern  varies  with  UT. 

•  Using  more  Inputs  and  being  subject  to  more  constraints  than 
ever  before,  a  theoretical  model  of  the  high-latitude  ionosphere 
was  able  to  successfully  reproduce  auroral-zone  observations 
from  two  widely  spaced  longitudes. 

•  Theoretical  calculations  of  the  6300-A  intensity  from  atomic 
oxygen  showed  that  the  hot  electrons  found  near  the  afternoon 
convection  reversal  are  responsible  for  a  much  larger  proportion 
of  the  emission  than  previously  considered. 

•  Although  convection  has  the  dominant  role  in  determining  the 
large-scale  nighttime  meridional  neutral  wind  at  F-region  alti¬ 
tudes,  local  effects  such  as  particle  and  Joule  heating  can 
create  significant  perturbations.  At  midnight  in  the  polar  cap 
a  new  phenomenon  was  observed,  an  abatement  or  sometimes  even  a 
reversal  of  the  southward  meridional  wind. 

•  The  importance  of  ion  drag  was  further  demonstrated  in  statis¬ 
tical  analyses  that  showed  the  meridional  component  of  the 
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neutral  wind  to  be  significantly  correlated  with  the  By  compon¬ 
ent  of  the  IMF.  These  analyses  also  showed  that  the  neutral 
wind  varies  significantly  with  season. 

•  Several  phenomena  observed  in  the  early  morning  at 
Chatanika — such  as  elevated  ion  and  neutral  temperatures, 
significantly  perturbed  meridional  neutral  winds,  and  reduced 
plasma  densities — could  be  the  result  of  very  large  Joule 
heating  Just  poleward  of  the  radar  or  in  the  polar  cap. 

•  New  methods  were  devised  to  determine  the  neutral  exospheric 
temperature  from  high  latitude  radar  data;  the  daily  and  solar- 
cycle  variation  were  studied. 

The  MITHRAS  observations  were  obtained  specifically  to  take  advan¬ 
tage  of  a  unique  radar  configuration  that  existed  for  a  short  period  in 
the  auroral  zone.  The  results,  to  date,  attest  to  the  importance  of  the 
data  set  as  well  as  to  the  effectiveness  of  the  research  approach. 
Additional  results  are  still  to  be  obtained  from  studies  begun  later  or 
based  on  the  results  of  earlier  ones. 

In  addition  to  what  can  be  learned  from  this  unique  data  set,  new 
possibilities  arise  from  the  current  configuration  of  the  radars. 
Sondrestrom,  EISCAT,  and  Millstone  Hill  are  now  located  so  as  to  provide 
denser  sampling  in  a  90*  segment  of  longitude  and  to  cover  a  wider 
latitude  region  from  the  polar  cap  to  Just  equatorward  of  the  auroral 
region.  This  configuration  offers  possibilities  for  more  detailed 
studies  of  many  aspects  of  the  interactions  among  the  magnetosphere, 
ionosphere,  and  neutral  atmosphere. 
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(May  1981  through  June  1982) 
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electric  field  response  to  substorms.  Radio  Scl. .  18.  981-987,  1983. 
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de  la  Beaujardiere,  0.,  V.  B.  Wickwar,  and  J.  D.  Kelly,  Observations  of 

polar-cap  ionospheric  signatures  of  solar  wind/magnetosphere  coupling. 
Proceedings  of  ARCAD  Workshop.  Toulouse.  M.  France,  1985. 

de  la  Beaujardiere,  0.,  V.  B.  Wickwar,  J.  D.  Kelly,  and  J.  H.  King,  IMF-By 
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116 


s;»  V 


.  c'i.,  f 


v  '  ■>  V  V->Vs>v*-  %•’,  .'.V>  '•*  V 


de  la  Beaujardlere,  O. ,  V.  B.  Wlckwar,  and  J.  H.  King,  Sondrestrom  radar 
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Wlckwar,  V.  B. ,  and  W.  Kofman,  Dayside  red  auroras  at  very  high  latitudes: 
The  importance  of  thermal  excitation,  Geophvs.  Res.  Lett..  2,  923-926, 
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118 
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